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The  energy  spectrum  of  a  5  curie  239Pu-Be  source  is 
determined  by  constructing,  calibrating  and  using  a  two 
detector  neutron  t ime-of -f 1 i ght  spectrometer.  This 
spectrometer  measures  the  time  it  takes  for  a  neutron  to 
travel  from  detector  A  to  point  B  following  an  elastic  scatter 
from  a  hydrogen  atom.  If  the  distance  of  travel  is  known, 
then  the  energy  of  the  neutron  is  then  calculated  using 
kinematics.  The  different  aspects  of  the  construction, 
circuit  operation,  and  experimental  geometry  are  explained. 

The  purpose  of  the  construction  of  the  spectrometer  is  for  use 
in  the  Nuclear  Instrumentation  Lab  at  the  Air  Force  Institute 
of  Technology.  The  time  spectrum  obtained  from  the 
spectrometer  is  unfolded  into  the  neutron  energy  spectrum.  A 
comparison  is  made  to  the  spectrum  published  by  Arun  Kumar  and 
P.  Nagarajan  and  others.  The  error  analysis  of  the  resolved 
energy  spectrum  is  calculated  using  a  Monte  Carlo  simulation 
routine  which  takes  1000  random  points  for  the  generation  of 
the  neutron  in  the  source,  and  for  the  interactions  in  the 
sensitive  volumes  of  both  detectors.  A  dual  parameter 


analyzer  is  used  to  look  at  the  coincidence  response  of  the 
recoil  protons  from  the  first  detector  to  the  time  spectrum 
from  the  experiment. 


Construction,  Calibration,  and  Use  of  a 
Neutron  T ime-of -F 1 i ght 
Spectrometer 


Since  Chadwick  first  discovered  the  neutron  in  1932, 
scientists  and  physicists  have  been  trying  to  determine  the 
energy  spectra  of  different  neutron  sources.  Several  types  of 
spectrometers  have  been  developed  through  the  years;  however, 
only  a  few  have  an  energy  range  wide  enough  to  be  practical. 
The  t ime-of -f 1 i ght  system,  which  measures  the  time  it  takes 
for  the  neutron  to  travel  from  point  A  to  point  B,  has  emerged 
as  the  preferential  choice  for  determining  the  energy 
spectrum.  This  system  should  provide  the  energy  spectrum  from 
any  source,  such  as  that  emitted  by  spontaneous  fission,  by  an 
(a,  n)  reaction,  or  by  any  nuclear  reaction.  If  the  time  for 
the  travel  of  a  fixed  distance  is  known,  then  the  kinetic 
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energy  of  the  neutron  can  be  calculated  through  kinematics. 
This  report  considers  the  different  aspects  of  the 
construction  of  this  type  of  spectrometer ,  the  calibration  and 
timing  of  the  device,  and  its  use  to  resolve  the  energy 
spectrum  of  a  5  curie  239Pu-Be  neutron  source. 


The  measurement  of  neutron  energies  with  high 
resolution  and  somewhat  reasonable  efficiency  is  difficult 
for  many  reasons.  The  biggest  factor  contributing  to  this 
difficulty  is  the  lack  of  a  net  electric  charge  on  the 
neutron,  which  eliminates  the  use  of  all  forms  of 
electromagnetic  and  direct  ionization  spectrometer s . 
Several  methods  have  been  developed  over  the  years  to 
measure  the  spectrum  of  different  neutron  sources.  Among 
these  are:  (1)  the  measurement  of  proton  recoil  energies; 
(1:126)  (2)  the  measurement  of  the  energy  release  in  a 
known  nuclear  reaction;  (2:380)  (3)  diffraction  of 
low-energy  neutrons  in  crystals;  (3:540)  (4) 

time-of-f 1 ight  measurements;  (3:541)  (5)  threshold 

detectors;  (4:708)  and  (6)  measurement  of  attenuation  and 
moderation  of  neutrons  in  various  thickness  of  hydrogenous 
material,  otherwise  known  as  the  use  of  Bonner  Spheres. 
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(6:95)  Each  of  these  methods  has  advantages  and 
disadvantages  and  each  class  is  useful  over  different 
ranges  of  neutron  energies. 

At  neutron  energies  above  100  keV,  it  is  possible  to 
detect  neutrons  with  reasonable  efficiency  by  detection  of 
recoiling  nuclei.  Recoil  protons,  which  are  generated 
from  hydrogenous  material,  are  measured  by  both  ionization 
and  range  measurements.  If  the  proton  recoil  energies  are 
measured  irrespective  of  direction  (4n),  then  the 
spectrum  of  the  protons  must  be  "unfolded"  by  use  of  a 
computer  program.  However,  if  the  proton  recoil  is 
measured  only  for  the  scatter  of  neutrons  in  a  fixed 
direction  relative  to  the  direction  of  the  incident 
neutron,  then  the  neutron  energy  is  directly  calculable 
from  the  energy  of  the  recoiling  protons.  The  main 
disadvantage  of  this  type  spectrometer  is  associated  with 
the  detection  of  neutrons  below  0.5  MeV.  This  limitation 
rises  from  the  minimum  track  length  for  detection  of 
proton  recoil  is  about  3  microns.  (3:540) 

Soon  after  the  discovery  of  the  neutron,  the  low 
efficiency  and  difficulty  in  detecting  neutrons  with 
reasonable  resolution  resulted  in  the  development  of  a 
number  of  spectrometer s  based  on  a  neutron- i nduced 
reaction  of  the  type  n+A^B+C  +  Q  .  For  these  reactions 
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the  product  particles  are  accompanied  by  an  energy 
release,  Q.  This  energy,  along  with  the  kinetic  energy  of 
the  neutron,  appears  as  the  kinetic  energy  of  the  product 
particles,  provided  they  are  not  in  an  excited  state. 

Even  though  the  individual  energies  of  particles  B  and  C 
will  vary  with  their  angles  of  emission,  the  total  energy 
is  always  the  sum  of  the  energy,  Q,  and  the  kinetic  energy 
of  the  incident  neutron.  The  neutron  energy  is  therefore 
obtained  from  a  measurement  of  the  total  energy  released. 

For  this  type  of  spectrometer  to  be  accurate,  the 
reaction  must  satisfy  several  requirements.  The  cross 
section  of  the  material.  A,  must  be  sufficiently  large  and 
a  smooth  function  of  the  neutron  energy.  The  products  of 
the  reaction  must  not  vary  in  the  decay  of  the  excited 
atom  which  is  generated  when  the  neutron  reacts  with  the 
nucleus.  Additionally,  if  the  low-energy  range  is  to  be 
considered,  the  value  of  the  energy  released,  Q,  must  be 
positive  and  the  material,  A,  must  be  of  a  form  that  is 
suitable  for  use  in,  or  with,  an  energy  sensitive 
detector.  Currently,  the  only  materials  which  come  close 
to  this  are  the  light  nuclei  3He,  6Li  and  10B. 

The  cross  section  for  absorption  of  neutrons  of  3He  is 
about  5327  barns  for  2200  m/sec  neutrons  and  varies  as  1/v 
for  neutrons  up  to  approximately  6  MeV.  (5:941)  The 


problem  associated  with  this  detector  system  is  the 
competing  reactions  which  reduce  the  number  of  (n,p) 
reactions.  The  most  significant  of  these  is  the  elastic 
scatter  cross  section  for  the  nucleus  of  the  helium  atom. 
This  scattering  cross  section  is  always  higher  than  the 
absorption  cross  section,  and  it  becomes  more  prominent  at 
higher  neutron  energies.  For  example,  the  cross  sections 
for  absorption  and  scatter  are  about  equal  at  a  neutron 
energy  of  150  keV,  but  the  scattering  cross  section  is 
about  three  times  more  probable  than  absorption  at  2  MeV . 
Additionally,  a  competing  (n,d)  reaction  on  3He  is 
possible  at  neutron  energies  exceeding  4.3  MeV,  but  this 
cross  section  is  low  for  neutron  energies  below  10  MeV. 
(4:527).  Therefore,  the  upper  limit  on  the  energy  of  the 
neutrons  which  it  can  effectively  detect  is  limited. 

A  widely  used  spectrometer  for  slow  neutrons  is  the 
BF3  proportional  tube  system.  This  device  uses  boron 
trifluoride  as  the  proport ional -counter  gas  and  10B  as  the 
target  for  converting  the  slow  neutrons  into  secondary 
particles.  A  number  of  other  boron-containing  gasses  have 
been  evaluated,  but  BF3  is  the  universal  choice  because  of 
its  superior  properties  as  a  proportional  gas,  as  well  as 
its  high  concentration  of  10B.  In  most  cases,  the  gas 
used  for  the  proportional  gas  is  highly  enriched  with  ,0B, 
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so  that  the  efficiency  is  some  five  times  greater  than  if 
the  gas  contained  natura 1 1 y-occurr i ng  boron.  This  type 
detector  gives  a  pulse  height  spectrum  for  the  bombarding 
neutrons  which  is  related  to  the  energy  spectrum  of  the 
neutrons . 

A  disadvantage  arises  with  the  use  of  the 
proport ional -tube  system  because  the  output  pulse  from  a 
BF3  tube  originates  with  reaction  products  created  with  a 
random  location  and  direction.  This  causes  typical  pulses 
with  rise  times  that  vary  by  as  much  as  3-5  microseconds 
for  tubes  of  average  size.  (4:493)  Another  disadvantage 
is  that  the  point  of  interaction  cannot  be  precisely 
defined  except  as  occurring  somewhere  in  the  volume  of  the 
gas.  With  a  typical  tube  being  10-20  centimeters  long  to 
provide  reasonable  interaction  efficiency,  path-length 
uncertainties  can  be  quite  large.  Therefore,  this  type  of 
system  is  rejected  for  the  application  presented  here. 

Another  widely  used  material  for  the  detection  of 
neutrons  is  6Li.  The  neutron  absorption  cross  section  for 
a  2200  m/sec  neutron  is  about  940  barns  and  varies  as  1/v 
for  neutrons  up  to  about  1  MeV  with  a  pronounced  resonance 
at  250  keV.  (5:941)  With  this  cross  section,  the  use  of 
6Li  becomes  a  very  capable  detector  for  low  energies; 
however,  for  fast  neutrons  over  1  MeV  the  material  becomes 
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ineffective.  A  competing  reaction,  6Li  (n,n’d)  4He,  has  a 
smaller  Q  value  and  becomes  the  dominant  neutron- i nduced 
reaction  above  about  2.5  MeV .  Because  this  reaction  leads 
to  three  products,  one  of  which  is  a  neutron  that  normally 
escapes,  a  continuum  of  deposited  energy  should  be 
expected  even  for  mono-energetic  incident  neutrons. 
Therefore,  even  though  it  will  contribute  some  neutron 
pulses,  this  additional  reaction  is  generally  an 
undesirable  part  of  the  response  of  any  detector  that 
attempts  to  measure  the  incident  neutron  energy.  (4:524) 
Crystal  Diffraction  neutron  spectrometers  with  high 
resolution  make  use  of  the  reflection  of  an  intense 
wel 1 -col  1 imated  beam  of  neutrons  at  the  Bragg  angles 
associated  with  the  atomic  planes  of  a  crystal.  A  highly 
collimated  beam  of  neutrons  from  a  nuclear  reactor  is 
allowed  to  fall  on  a  large  crystal  at  an  angle  e  relative 
to  the  crystal  planes.  Some  of  the  neutrons  which  have  a 
wave  length  A.  such  that  they  satisfy  the  Bragg 
relationship  are  reflected  through  an  exit  collimator  into 
a  neutron  detector  placed  at  the  angle  0  relative  to  the 
reflecting  planes.  This  relationship  for  the  detector  is: 


where  n  is  an  integer,  X  is  the  DeBroglie  wavelength  of 
the  neutron,  d  is  the  spacing  between  crystal  planes  and  0 
is  the  angle  between  the  direction  of  the  incident  neutron 
beam  and  the  reflecting  plane.  The  resolving  power  of 
this  spectrometer  is  largely  governed  by  the  uncertainties 
in  defining  the  angle  0  .  Additionally,  there  is  an 

intrinsic  uncertainty  introduced  by  imperfections  in  the 
crystal  structure  itself,  and  there  is  an  ambiguity 
induced  because  n  can  have  values  of  1,2,3,  etc.  Thus  if 
X2  is  equal  to  0.5X,  and  nT  -  2  and  n,  =1  then  0,  =  02 
and  the  accuracy  of  the  resolution  is  gone.  Crystal 
spectrometers  have  been  used,  however,  with  great  success 
up  to  energies  of  about  10  eV.  Their  performance 
deteriorates  above  this  due  to  a  loss  in  the  counting  rate 
as  the  neutron  energy  increases.  (3:540-541) 

None  of  the  spectrometer s  mentioned  thus  far  cover  all 
energy  ranges.  There  is,  however,  one  type  which  will 
cover  all  energy  ranges  of  interest,  the  t ime-of -f 1 i ght 
spectrometer.  Time-of-f 1 i ght  spectrometer s  measure  the 
flight  time  between  two  fixed  points.  There  are  several 
different  types  of  this  spectrometer .  The  first  of  these 
is  the  pulse  t ime-of-f 1 i ght  spectrometers,  in  which  the 
neutrons  are  pulsed  out  of  the  source  and  the  time  it 
takes  for  the  different  energy  neutrons  to  travel  a  fixed 
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distance,  l  ,  is  measured.  This  time,  A t  ,  is  collected 
and  the  time  spectrum  converted  into  the  energy  spectrum 
of  the  neutron  source.  Not  having  a  pulse  accelerator  at 
all  facilities  led  into  the  development  of  a  different 
type  of  t ime-of -f 1 i ght  spectrometer . 

This  type  of  spectrometer  uses  the  scatter  of  a 
neutron  from  a  hydrogenous  material  into  a  second  detector 
a  fixed  distance  away  from  the  first  detector.  With  this 
system,  the  time  of  the  flight  of  a  scattered  neutron  from 
the  first  detector  to  the  second  detector  can  be  measured. 
Figure  (1)  is  a  typical  set  up  for  the  two  detector 
system.  With  the  system  configured  this  way  the  neutron 
could  travel  the  distance  from  the  source  to  the  first 
detector,  scatter  into  the  solid  angle  of  the  second 
detector  and  generate  both  start  and  stop  times  for  a 
device  which  could  measure  this  flight  time.  The  use  of  a 
t ime-to-ampl i tude  converter  does  this  job.  By  setting  the 
time  range  of  the  converter  to  a  maximum  time  in  which  the 
slowest  neutron  of  interest  (0.5  MeV)  would  take  to 
traverse  the  distance  between  the  detectors  signals  of 
different  pulse  heights  would  be  generated  for  the  various 
flight  times.  These  signals  could  be  fed  into  a 
Multi-Channel  Analyzer  and  the  time  spectrum  generated 
over  a  period  of  time.  Care  must  be  taken  when  the  time 
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spectrum  is  unfolded  into  the  energy  spectrum  of  the 
neutron  source  in  order  to  properly  calculate  the  incident 
neutron  energy.  For  scattered  neutrons  which  are  not 
relativistic,  this  t ime-of - f 1 i ght ,  At  ,  is  given  in 
nanoseconds  by  the  relation 


At 


72.31 


(2) 


where  l  is  the  length  of  the  flight  path  in  meters  and  E  is 
the  energy  of  the  neutron  in  MeV.  By  determining  the 
values  of  At  over  a  fixed  distance  l,  the  energy  spectrum 
for  any  neutron  source  can  be  determined.  Should  the 
neutrons  approach  the  rest  mass  of  the  neutron  (939.566 
MeV),  equation  2  is  no  longer  valid  and  the  kinetic  energy 
of  the  neutron,  E „,  can  be  determined  through  the  use  of 
the  familiar  result  of  special  relativity:  (1:541) 
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where  E0  is  the  rest  mass  of  the  neutron  and  c  is  the  speed 


of  light.  Substituting  v-l/At  yields 


If  the  units  of  energy  are  MeV,  and  those  of  length 
and  time  are  meters  and  nanoseconds,  respectively,  then 
the  time  it  takes  in  nanoseconds  for  the  neutron  to  travel 
distance  l,  is 


A  t 


3.3356410 1 


fo  yV 


(5) 


Since  the  time  of  travel  can  be  measured  with  a  proper 
geometry  set-up  of  equipment  and  the  proper  electronics, 
the  energy  of  the  scattered  neutron  £ ,  can  be  calculated. 
Therefore,  the  energy  limitations  which  existed  on  the 
other  types  of  spectrometer s  do  not  exist  with  the 
t ime-of -f 1 i ght  system  which  is  it's  biggest  advantage. 
Another  advantage  of  the  t ime-of -f 1 i ght  system  is  the 
simplicity  of  the  electronics  involved.  The  use  of  a 
t ime-to-amp 1 i tude  converter  to  measure  the  flight  time  is 
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NEUTRON  TROM  SOURCE  ENTERS  DETECTOR  AT  POINT  A 
SCATTERS  THROUGH  ANGLE 

TRAVELS  DISTANCE  I  AND  ENTERS  DETECTOR  AT  POINT  B 


Figure  1.  Typical  Two  Detector 
Time-of-F 1 i ght  Spectrometer  Geometry 


vital  to  the  ease  in  which  the  time  spectrum  can  be 
generated  for  the  two  detector  system.  The  greatest 
disadvantage  of  the  system  is  its  inefficiency.  This 
inefficiency  comes  about,  for  small  neutron  sources,  from 
the  number  of  scattered  neutrons  which  are  incident  on  the 
second  detector,  resulting  in  long  counting  times  in  which 
the  linearity  of  the  counting  system  may  fail. 

The  primary  concern  for  the  measurement  of  neutron 
spectra  is  the  energy  range  that  the  method  will  cover. 
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Figure  (2)  is  a  graphical  representation  of  the  energy 
ranges  that  can  be  effectively  analyzed  by  the  different 
detection  methods.  (6:96) 


0.001  0.01  0.1  1  10  100 

Neutron  Energy  (MeV) 

Figure  2.  Energy  Ranges  for  Neutron  Spectrometers 

As  depicted  in  figure  (2),  the  best  overall  energy 
spread  is  obtained  with  either  a  t ime-of-f 1 i ght 
spectrometer  or  with  a  Bonner  sphere  system.  The 
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limitations  of  the  energy  range  of  the  other  detectors 
does  not  lend  itself  to  this  study  of  the  energy  spectrum 
from  a  neutron  source  such  as  239Pu-Be,  which  is  the 
source  of  neutrons  for  analysis  in  this  report. 

Scope 

The  scope  of  this  project  deals  with  the  different 
aspects  of  the  t ime-of-f 1 i ght  spectrometer  setup, 
calibration,  and  use.  The  main  objective  of  this  study  is 
to  construct,  test,  calibrate  and  use  a  neutron  time  of 
flight  detector.  This  is  accomplished  through  the  use  of 
experimentation  along  with  the  unfolding  of  the  time 
spectrum  into  the  energy  of  the  incident  neutrons.  A  5 
curie  239Pu-Be  source  will  be  used  to  test  the  detectors 
and  prove  out  the  electronics.  The  energy  range  of  this 
source  is  from  0.5  MeV  to  10.5  MeV.  Therefore  classical 
kinematics  will  be  used  for  the  resolution  of  the  energy 
spectrum  from  the  time  spectrum. 

Assumptions 

In  order  to  complete  the  study,  the  following 
assumptions  are  made: 
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(1)  The  response  of  all  electronic  equipment  is 
1 i near . 

(2)  During  the  t ime-of - f 1 i ght  measurement,  the  fast 
neutrons  will  not  overtake  slow  neutrons  causing  errors  to 
be  generated  in  the  time  spectrum. 

(4)  All  measurements  of  length  in  the  lab  will  be 
accurate  to  =*=  2%. 

(5)  Multiple  scatters  in  the  sensitive  volumes  of  the 
detectors  can  be  neglected  which  would  degrade  the  energy 
of  the  incident  neutron  yielding  a  softer  energy  spectrum. 

(6)  The  neutrons  will  travel  in  a  straight  path 
between  the  two  detectors. 

(7)  Even  though  the  plastic  scintillators  have  almost 
as  many  carbon  atoms  as  hydrogen  atoms,  the  scattering  of 
the  neutrons  from  these  atoms  will  not  be  considered. 

Sequence  of  Presentation 

The  theory  of  the  unfolding  of  the  time  spectrum  into 
the  energy  spectrum  of  the  incident  neutrons  is  discussed, 
along  with  the  factors  that  deal  with  the  detection 
efficiency.  Following  this,  the  equipment  selected  for 
the  detection  system  is  described  along  with  the  operation 
of  those  items  which  may  not  be  familiar.  The  procedure 
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for  the  experimentation  is  explained  in  detail  and  the 
results,  along  with  the  error  analysis,  provided. 
Finally,  the  conclusions  are  reached  concerning  the 
detection  system  and  the  analysis  of  the  data. 
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II.  Experimental  Considerations 


To  successfully  undertake  this  assignment  a 
spectrometer  is  constructed  to  process  the  signals  from 
the  anodes  of  two  plastic  scintillator  detectors  which 
correspond  to  the  time  difference  between  a  start  signal 
and  a  stop  signal  generated  by  the  same  neutron.  Constant 
fraction  d  i  scr  imi  nator  s  ,  («.FD),  are  used  to  reduce  the 
signal  jitter  from  the  anodes  and  to  eliminate  the 
dependence  of  the  amplitude  o,  the  signal.  The  start  and 
stop  signals  are  fed  into  a  t ime-to-amp 1 i tude  converter 
which  will  generate  a  0  to  10  volt  logic  signal,  depending 
on  the  time  between  start  and  stop  signals  and  the  time 
range  selected  for  analysis.  A  time  spectrum  is  generated 
and  this  spectrum  unfolded  to  produce  the  energy  spectrum 
of  the  source.  This  system  is  set  up  initially  in  a  high 
scatter  environment  with  the  use  of  a  neutron  moderator  to 
form  a  collimated  beam  of  neutrons  from  the  source.  The 
geometry  of  the  detector  set  up  for  experiments  1  and  2  is 
shown  in  figure  (3).  Experiment  1  will  have  the  detectors 
1.5  meters  apart  and  experiment  2  will  place  the  detectors 
0.65  meters  apart. 
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Figure  3.  Experiment  Geometry 

i 

The  copper  shield  reduces  the  number  of  first  flight 
incident  neutrons  from  the  source  on  the  second  detector. 

,  The  shadow  shield  helps  attenuate  the  number  of  second  and 

third  scatter  neutrons  from  the  surroundings  of  the 
detectors.  A  discussion  on  the  construction  of  the  shadow 

i 
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shield  is  given  in  a  later  section  of  the  document. 
Following  the  experiments  within  the  high  scatter 
environment,  the  detectors  and  electronics  are  moved  into 
a  low  scatter  environment  and  the  experiments  rerun. 

After  a  detailed  analysis  of  the  data  obtained  from  the 
first  four  experiments,  the  plastic  scintillators  are 
replaced  with  ones  that  have  twice  the  volume  of  the  first 
ones  and  the  experiments  conducted  again.  With  this 
complete,  an  analysis  is  done  using  the  dual  parameter 
feature  of  the  analyzer,  by  placing  the  recoil  proton 
spectra  on  one  axis,  the  t ime-of-f 1 i ght  information  on  the 
other  axis,  and  determining  the  relationship  between  the 
proton  recoil  energy  and  the  t ime-of-f 1 i ght  information. 

Theory 

The  detection  of  high  energy  neutrons  has  challenged 
scientists  and  physicists  since  the  early  1930’s.  The 
most  common  method  of  determining  the  presence  of  these 
neutrons  is  based  on  the  scattering  of  the  neutrons  by 
light  nuclei.  Elastic  scattering  transfers  some  portion 
of  the  energy  into  the  target  nucleus,  which  recoils  and 
acts  much  like  a  proton  or  alpha  particle  as  it  loses  its 
energy  in  the  detector  medium.  Because  of  its  relative 
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size  to  that  of  the  neutron  and  its  behavior  as  a  recoil 


proton,  hydrogen  is  the  preferred  choice.  To  maximize  the 
number  of  hydrogen  molecules  per  cubic  centimeter,  plastic 
scintillators  were  selected  to  generate  the  recoil  proton 
and  subsequent  signals  from  the  interaction  of  the 
neutrons . 

To  begin  the  signal  generation  in  the  sensitive  volume 
of  the  scintillators,  the  target  nuclei  are  for  all 
practical  purposes  at  rest  prior  to  the  collision.  During 
a  single  hydrogen  collision,  the  amount  of  energy  transfer 
to  the  nucleus  can  be  from  0  to  the  full  energy  of  the 
neutron,  which  results  in  a  signal  proportional  to  the 
energy  deposited.  To  examine  this  better,  the  kinematics 
of  neutron  elastic  scattering  must  be  investigated. 

A  center  of  mass  system  with  non-re  1  at i v i st i c  neutrons 
(En  <<  939  MeV)  is  assumed.  By  conserving  momentum  and 
energy  we  are  able  to  determine  the  energy  of  the 
scattered  neutron  from  the  inelastic  collision  with  a 
proton  by  use  of  the  equation 


2/1 


(1  +/1): 


;(  ]  -  CO S0)/f 


N 


(6) 
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where  A  is  the  mass  of  the  target  nucleus,  En  is  the 
energy  of  the  incoming  neutron,  Es  is  the  energy  of  the 
scattered  neutron,  and  6  is  the  scattering  angle  of  the 
neutron  in  the  center  of  mass  system. 

To  convert  this  equation  into  a  laboratory  system 
where  the  protons  are  initially  at  rest,  we  use  the 
transf ormat ion 


COS0  = 


1  -  COS0 


2 


(7) 


where  e  is  the  scattering  angle  of  the  recoil  nucleus. 

Combining  the  last  two  equations  we  have 

4A  ?  (8) 

-~-2Cc°s70)Fn 

( 1  +  Ay 

which  for  neutrons  which  scatter  from  hydrogen  atoms 
reduces  to 
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Fs  -  (cos zQ)En 


(9) 


Scatters  of  neutrons  from  carbon  atoms,  within  the 
sensitive  volume  of  the  detector,  will  impart  only  14%  of 
there  energy  into  the  scintillator  which  will  make  the 
energy  spectrum  below  1.5  MeV  suspect  to  the  final 
resolution.  This  results  from  calculation  of  equation  (8) 
using  the  value  of  12  for  the  atomic  number  of  carbon  and 
assuming  a  scatter  angle  of  45° .  Therefore,  the 
assumption  that  the  carbon  atoms  will  not  contribute  to 
the  final  resolution  is  somewhat  valid. 


To  determine  the  value  of  Es  ,  classical  kinematics 
along  with  the  application  of  distribution  function 
manipulation  is  used.  To  do  this  we  start  with  the 
equation  of  kinetic  energy  for  the  neutron 


1  2 
*  - mv 
5  2 


(10) 
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If  the  neutron  travels  from  the  first  detector  to  the 
second  detector  a  known  distance  l  in  time  At,  then  the 
velocity  of  the  neutron  is 


v  “ 


At 


(11) 


By  differentiating  the  equation  with  respect  to  A?  we 

get 


dv  _  l  (12) 

dCAn~~At* 

or  that  dv  =  -l/Atz  dt.  To  eliminate  dependance  on  the 
number  of  particles  the  spectrum  of  the  time,  velocity, 
and  energy  will  be  reduced  to  a  normalized  spectrum.  To 
accomplish  this,  the  time  spectrum,  which  is  measured  in 
the  laboratory,  is  divided  by  the  total  integral  of  the 
spectrum  resulting  in  S( Af)  .  This  relationship  is  needed 
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for  the  next  step  in  the  derivation,  which  is  to  go  from 
the  time  distribution  to  the  velocity  distribution  and 
finally  to  the  energy  distribution  through  the  use  of 


S(t;)  =  S(AO 


d(AQ 

dv 


(13) 


which  leads  to 


S(p)  -  S(At) 


A  t? 

~T 


(14) 


Correlating  the  energy  to  the  velocity  gives 


dv 

S(/Ts)  =  S(u)  — 
d  E 


(15) 


Returning  to  equation  (9)  we  find  the  alliance  between  the 
energy  of  the  particle  to  the  velocity  d.E,/dv  of  the 
particle  to  be 
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mu 


(16) 


dEs 

dv 


and 


dv  1 

dEs  *  ffi 

dv 


I  _  At 
mu  mZ 


(17) 


substituting  the  previous  values  for  v,  dv ,  and  realizing 
that  the  value  of  dv  is  the  absolute  value  of  the 
relationship  in  equation  (12),  we  have 


S(£s)«S(t/) 
S(u)-S(A() 
S(fs)-S(At) 
5(£s) 


S(At) 


du 


d£. 


d&<  du 


dv 
dE~s 
d(AQ 
dp 

dAf  dp 
dv  dE s 
£(AQ 


(18) 


Af 3  (19) 

S(£s)-  —  S(At) 
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or  that  the  distribution  of  particles  of  the  scattered 
energy  is  equal  to  the  constants  in  equation  (19)  times 
the  value  of  the  distribution  in  time  t.  The  energy  of 
the  source  neutron  is  found  by  normalizing  the  energy 
spectrum  from  the  relationship  in  equation  (9). 

S(Fs)-S(£w)Cos2G  (20) 

which  converts  to  the  total  transformation  of  the  time 
spectrum  into  the  incident  neutron  energy  spectrum  by 
using 


S(£„)-S(  AO 


At' 


(21) 


ml2  Cos?G 

The  intrinsic  detection  efficiency  of  any  detector 
based  on  recoil  nuclei  can  be  calculated  from  the 
scattering  cross  section  o,  .  For  the  Bicron 
scintillators  we  have  both  carbon  and  hydrogen  atoms,  and 
the  competing  effects  of  the  carbon  scattering  must  be 
taken  into  account.  Utilizing  the  relationship  in 
Equation  (8),  for  the  carbon  atoms  with  an  atomic  number 
of  12,  and  a  scatter  angle  of  45°  ,  results  in  a  transfer 
of  14  %  of  the  energy  of  the  incident  particle  into  the 
carbon  atom.  To  ensure  the  data  from  the  recoil  of  the 
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neutron  from  the  hydrogen  atoms  is  not  contaminated,  the 
Lower  Level  Discriminator,  (LLD),  must  be  set  to  eliminate 
the  carbon  scatters.  The  intrinsic  efficiency  for  an 
axial  beam  of  neutrons,  neglecting  multiple  scatters 
within  the  sensitive  volume  of  the  detector,  is  given  by 


N„  a„  (22) 

6-  -(l-exp(-(A/Ha„-WVcoc)d)) 

where  the  subscripts  H  and  C  refer  to  the  separate 
hydrogen  and  carbon  quantities,  and  d  is  the  path  length 
of  the  neutron  travel.  (5:534) 

Equation  (22)  can  be  explained  in  the  following 
manner.  The  first  term  of  the  equation  is  the  probability 
of  the  incident  neutron  coming  in  contact  with  a  hydrogen 
atom  vice  the  carbon  atom.  The  second  term  is  the 
probability  that  a  neutron  which  travels  a  distance,  d  , 
in  the  sensitive  volume  medium  will  have  an  interaction 
with  the  medium.  Therefore,  the  probability  of 
interaction  with  a  hydrogen  atom  times  the  probability  of 
interaction  within  the  medium  converts  to  the  intrinsic 
efficiency  for  the  detector. 
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An  empirical  fit  to  the  hydrogen  scattering  cross 


section  given  by  Marion  and  Young  is 


<*s  (*„) 


4-83 

-  0.578  barns 

{En 


(23) 


where  the  energy  of  the  incident  neutron  Fn  is  in  MeV 
(5:50/lN.  The  cross  section  for  the  carbon  atoms  in  the 
rr  on  of  interest  has  many  resonance  peaks  in  the  area  of 
the  energy  spectrum  and  would  have  to  be  individually 
calculated.  For  this  study  the  average  cross  section  of 
the  carbon  atoms  is  placed  at  2  barns.  With  this  cross 
section,  the  probability  of  entering  into  a  collision  with 
the  carbon  atoms  is  greater  than  for  the  hydrogen  atoms. 
This  leads  to  a  lower  intrinsic  efficiency  of  the  system 
and  a  source  of  possible  erroneous  signal  generation  in 
the  detectors.  Most  of  this  is  overcome  by  the  adjustment 
of  the  Lower  Level  Di scr imi nator .  The  efficiencies  of  the 
detectors  will  be  calculated  and  reported  in  the  results 
section  of  this  report. 
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Equipmenj; 


The  detectors  and  electronics  used  for  the  detection 
system  were  selected  based  on  the  response  characteristics 
of  the  different  modules.  A  detailed  list  of  equipment 
used  in  the  detection  system  is  listed  in  Appendix  A.  The 
detection  circuit  is  shown  in  figure  3. 

Several  factors  went  into  the  selection  of  the 
different  components.  The  desirable  features  of  a  neutron 
detector  used  in  a  t ime-of - f 1 i ght  spectrometer  are: 

(3:559) 

1)  high  efficiency  which  is  a  slowly  varying  function 
of  energy, 

2)  low  efficiency  for  detecting  background  radiation, 

3)  fast  time  response,  preferably  in  the  nanosecond 
range , 

4)  ease  of  manufacture  in  various  shapes  and  sizes, 

5)  simplicity  of  the  associated  electronics. 

Based  on  these  features,  the  use  of  an  organic 
scintillator  was  selected  for  the  detectors.  The  first 
scintillator  used  was  a  Bicron  422  plastic  scintillator 
which  is  made  of  Po 1 y v i ny 1  to  1 uene  and  has  a  decay  time  of 
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CFD  Gate  and 
Delay  Amp 


Figure  4.  Detector  System 

Abbreviations  used  in  figure  (3)  include:  PA  =  Preamp, 
LA  =  Linear  Amp,  CFD  =  Constant  Fraction  Di scr imi nators , 
TSCA  =  Timing  Single  Channel  Analyzer,  TAC  = 

T ime- to-Amp 1 i tude  Converter. 
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1.6  nanoseconds.  The  scintillator  was  a  cylinder  2  inches 
in  diameter  and  1  inch  high  and  had  a  carbon  to  hydrogen 
molecular  ratio  of  0.91.  This  scintillator  was  used  for 
experiments  1-4.  Following  experiment  4  the  scintillators 
in  both  detectors  were  exchanged  for  a  Bicron  404  plastic 
scintillator  coated  with  a  titanium  dioxide  on  the  outside 
of  the  cylinder  with  the  two  flat  faces  clear.  This 
scintillator  is  also  made  of  Po 1 y v i ny 1  to  1 uene ,  has  a  decay 
time  of  1.8  nanoseconds,  and  is  the  same  diameter  as  the 
previous  ones  except  the  height  of  the  scintillators  is  2 
inches  and  has  a  carbon  to  hydrogen  molecular  ratio  of 
0.91  also.  This  doubled  the  volume  of  the  scintillators, 
and  with  the  hydrogen  concentration  of  the  two 
scintillators  nearly  the  same  (5.21  X  1022  atoms/cm3  vs 
5.19  X  1022  atoms/cm3),  the  number  of  neutron  interactions 
were  expected  to  double  on  the  first  detector  and  increase 
by  a  factor  of  four  on  the  second  detector.  The  signals 
generated  from  the  scintillators  will  vary  in  size 
depending  of  the  amount  of  energy  deposited  and  the  angle 
of  scatter  from  the  first  and  second  detectors.  A  timing 
device  was  used  to  reduce  this  dependence. 

Because  the  base  routes  negative  high  voltage  to  the 
photo-cathode  of  the  photo-multiplier  tubes  and  the  anode 
is  grounded,  the  detector  assembly  consisting  of  the 
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scintillator,  photo-multiplier  tube  and  base  had  to  be 
modified.  To  prevent  a  voltage  drop  across  the  glass 
envelope  of  the  photo-multiplier  tube  and  to  minimize 
accidental  shock,  aluminum  foil  was  used  to  cover  the 
scintillator  and  a  portion  of  the  tube.  A  wire  was 
soldered  from  the  high  voltage  input  at  the  bottom  of  the 
base  and  attached  to  the  aluminum  foil  with  black 
electrical  tape.  This  modification  kept  the  photo-cathode 
and  the  scintillator  at  the  same  negative  potential.  To 
minimize  loading  of  the  power  supply  in  case  of  a  short,  a 
10  Mf2  resistor  was  added  in  series  between  the  high 
voltage  input  and  the  aluminum  foil. 

Black  electrical  tape  was  utilized  to  prevent  the 
scintillators  from  coming  off  and  to  prevent  light  leaks 
which  could  render  the  photo-tubes  useless.  Furthermore, 
an  aluminum  cover  and  end  cap  were  slipped  over  the  tube 
assembly  and  grounded  to  the  base  structure  by  screws  to 
protect  the  user. 

Constant  fraction  d i scr imi nator s  were  used  to  reduce 
the  jitter  of  the  output  of  the  anodes  and  to  eliminate 
the  amplitude  dependence  on  the  output  of  the  signal. 

This  technique  is  particularly  useful  with  detectors  in 
which  the  collection  time  varies  according  to  the  location 
of  the  initial  ionization  in  the  detector.  The  constant 
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fraction  technique  operates  as  follows:  An  incoming 
unipolar  signal  is  split  into  two  paths.  The  signal  is 
attenuated  in  one  path,  and  in  the  other  path  the  signal 
is  not  attenuated,  but  delayed.  The  attenuated  and 
delayed  signals  are  then  combined  in  a  differential 
amplifier  stage,  resulting  in  a  subtraction  process.  This 
signal  is  bipolar,  with  the  baseline  crossing  being  the 
locus  of  a  particular  point  on  the  leading-edge  of  the 
incoming  signal.  This  point  is  at  a  fraction  of  the 
amplitude  of  the  incoming  signal,  and  that  fraction  does 
not  vary  with  amplitude.  A  crossover  pickoff  circuit 
generates  a  negative  logic  pulse  which  is  time  correlated 
to  the  desired  leading  edge  fraction.  These  devices 
generate  the  timing  signals  for  the  t ime-to-ampl itude 
converter . 

A  t ime- to-amp 1 i tude  converter  was  selected  based  on 
its  time  scales  and  time  resolution.  Input  signals  are 
generated  by  the  detection  system  to  give  a  start  and  stop 
signal  for  the  device.  Based  on  the  time  between  start 
and  stop  signals,  the  range  selected,  and  the  multiplier 
setting,  a  logic  signal  from  0  to  10  volts  is  generated. 
The  Tennelec  model  selected  has  a  range  selection  of  25 
nanoseconds  to  3  milliseconds.  For  a  neutron  source  with 
an  energy  range  of  0.5  to  10.5  MeV,  these  limitations 


33 


place  a  spread  in  the  distance  between  the  two  detectors 
such  that  0.6  <  l  <  95.074  meters.  With  this  large  of  a 
spread  in  the  values  of  l  ,  the  two  detectors  can  be  placed 
in  an  infinite  number  of  positions.  However;  the  further 
apart  the  detectors  are,  the  less  efficient  the  system 
becomes.  The  longer  flight  paths  lead  to  a  fewer  number 
of  neutrons  which  are  incident  on  the  second  detector 
through  spherical  divergence,  however;  the  longer  flight 
paths  gives  a  better  definition  of  the  scatter  angle  6.  To 
optimize  the  trade-offs  and  to  improve  the  efficiency  of 
the  spectrometer ,  the  maximum  distance  for  consideration 
was  placed  at  1.5  meters.  Time  resolution  of  the 
converter  is  listed  as  0.01&  of  full  scale  plus  5 
picoseconds . 

C ir cm it  Operation 

To  generate  a  time  spectrum,  the  electronics  of  the 
system  must  be  able  to  select  those  neutrons  which  are 
incident  on  the  first  detector  and  then  scatter  into 
the  second  detector  and  interact.  Should  the 
electronics  be  able  to  do  this,  then  the  time  spectrum 
can  be  generated  through  the  use  of  a 
t ime-to-pu 1 se-he i ght  converter  (TTPHC) . 
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Start  and  stop  signals  for  the  t ime- to-amp 1 i tude 
converter  must  be  generated  to  obtain  a  time  spectrum. 
To  obtain  these  signals,  a  neutron  from  the  source 
traverses  the  air,  enters  the  sensitive  volume  of  the 
first  detector,  has  an  elastic  collision  with  the 
first  detector  and  the  recoil  proton  generates  a 
negative  signal  on  the  anode  in  the  first  detector. 

If  the  same  neutron  scatters  in  the  direction  of  the 
second  detector  and  has  an  elastic  collision  with  the 
hydrogen  molecule  in  the  sensitive  volume,  it  will 
generate  a  signal  at  the  anode  of  the  second 
photo-multiplier.  Because  of  the  many  orders  of 
magnitude  of  neutrons  which  are  incident  on  the  first 
detector  and  scatter  in  many  different  angles,  the 
signals  from  the  first  detector  do  not  always  result 
in  a  signal  generation  in  the  second  detector,  which 
could  overwork  the  TTPHC.  To  eliminate  overworking 
the  TTPHC,  the  signal  from  the  first  detector  is 
delayed  through  the  use  of  a  fixed  time  delay  which 
allows  the  t ime- to- amp  1 i tude  converter  to  look  only 
for  response  pulses  from  the  neutrons  which  make  it  to 
the  second,  or  scattered  detector.  The  delay  time 
must  be  greater  than  the  amount  of  time  it  would  take 
for  the  slowest  neutron  to  travel  the  distance  between 
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the  two  detectors.  These  signals  from  the  photo-tube 
anodes  are  fed  to  the  constant  fraction  d i scr imi nator s 
which  generate  a  negative  logic  pulse.  From  the  fixed 
delay,  the  signal  from  the  first  detector  is  sent  to  a 
gate  and  delay  circuit  to  correct  for  differences  in 
the  signal  generation  timing  caused  by  the  electronics 
and  then  to  the  stop  side  of  the  time  to  amplitude 
converter.  The  signal  from  the  second  detector 
constant  fraction  di scr imi nator  is  also  sent  to  a  gate 
and  delay  circuit  for  timing  purposes,  and  then  to  a 
fast  coincidence  circuit  awaiting  the  response  of  the 
Pulse  Shape  analyzer. 

Since  to  the  differential  cross  section  of  the 
hydrogen  atoms  is  nearly  constant  for  scattering  into 
any  angle  0  ,  the  neutrons  scatter  isotropically  from 
the  first  detector,  of  which  a  certain  percentage  are 
incident  on  the  second  detector.  The  percentage  of 
scattered  neutrons  that  are  incident  on  the  second 
detector  depends  principally  on  the  distance  between 
the  first  detector  and  the  second  detector.  For  those 
neutrons  which  interact  with  the  second  detector,  a 
pulse  similar  to  that  from  the  first  detector  is 
generated.  To  prevent  any  direct  flight  neutrons  from 
contaminating  the  data,  a  copper  bar,  with 
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measurements  of  3"  X  4"  X  24",  was  placed  between  the 
second  detector  and  the  source  (See  figure  2).  The 
copper  bar  allows  the  neutron  to  scatter,  lose  energy 
and  become  absorbed  through  an  (n,y)  reaction.  These 
photons  could  interact  with  the  sensitive  volume 
principally  through  Compton  scatter  to  produce  a 
signal  which  would  effect  the  neutron  data. 
Interference  from  these  photons  is  eliminated  by  using 
a  Pulse-Shape-Analyzer.  This  device  filters  out  all 
signals  which  have  a  fall  time  of  less  than  1 
microsecond  from  the  dynode  of  the  second  detector. 

The  photons  produces  signals  that  have  fall  times  less 
than  this  and  the  neutrons  produce  signals  that  are 
greater  than  this  value.  Timing  for  the  fast 
coincidence  circuit  is  provided  for  the  output  signal 
of  the  Pulse  Shape  Analyzer  through  a  gate  and  delay 
timing  module.  A  signal  which  makes  it  through  the 
Pulse  Shape  Analyzer  is  from  a  neutron,  not  a  photon, 
and  this  signal,  when  in  coincidence  with  the  anode 
signal,  will  generate  the  start  signal  for  the 
t ime-to-amp 1 i tude  converter.  Therefore,  the  second 
detector  generates  the  start  signal  after  pulse  shape 
analysis  and  the  first  detector  generates  the  stop 
signal . 
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To  generate  proportional  pulse  heights 
corresponding  to  the  different  t ime-of -f 1 i ghts ,  the 
t ime- to-amp 1 i tude  converter  time  range  is  set  equal  to 
the  delay  time  of  the  fixed  delay  amp  in  the  first 
detector.  Measurements  of  TOF  were  made  for 
separation  distances  between  the  two  detectors.  Since 
a  lower  bound  of  0.5  MeV  was  selected  for  the  range  of 
energies  to  be  measured,  the  longest  flight  times  for 
the  neutrons  to  travel  1.5  and  0.65  meters  are 
respectively  153  nanoseconds  and  100  nanoseconds. 
Therefore,  the  signals  from  the  first  detector  was 
delayed  by  200  nanoseconds  to  allow  the  signal  from 
the  second  detector  to  start  the  TTPHC. 

To  generate  only  an  8  volt  signal  for  the 
multi-channel  analyzer,  the  multiplication 
potentiometer  was  set  at  1.25.  This  ensured  that  the 
maximum  voltage  for  a  200  nanosecond  travel  would  be  8 
volts. 

Procedure 

A  total  of  seven  experiments  were  conducted  to  study 
various  response  of  the  neutron  spectrometer .  In  the 
first  two  experiments  t ime-of - f 1 i ght  spectrum  were 


measured  in  a  high-scatter  environment,  and  with  the  use 
of  the  neutron  moderator  to  form  a  collimated  beam  of 
neutrons.  This  environment  consisted  of  concrete  floor 
along  with  concrete  block  walls  to  scatter  the  neutrons. 
Detectors  were  set  up  as  illustrated  in  figure  2,  and  the 
detector  circuit  connected  as  depicted  in  figure  3.  A 
negative  high  voltage  of  -2200  volts  was  placed  on  the 
photo-tubes  and  allowed  to  stabilize  for  twenty-four 
hours.  The  detection  system  was  timed  and  calibrated 
using  the  procedure  described  in  Appendix  B. 

The  source  selected  for  analysis  is  a  239Pu-Be  source 
obtained  from  the  Monsanto  Research  Corporation  in 
Miamisburg,  Ohio  in  1962.  The  initial  activity  strength 
of  the  239Pu  was  5  curies  with  a  total  output  of  9.04  X 
10*  neutrons  per  second  when  first  purchased.  However, 
because  the  source  is  contaminated  with  a  small  amount  of 
24 1 Pu,  the  number  of  neutrons  emitted  has  increased. 

Today  the  number  of  neutrons  that  the  source  produces  is 
1.1  X  107  per  second.  This  increase  is  caused  by  the 
decay  of  24 'Pu  to  2 4 1  Am  which  then  decays  by  the  emission 
of  an  a  particle.  With  the  increased  number  of  a  particles 
to  react  with  the  Beryllium,  the  number  of  neutrons 


Emission  of  neutrons  from  the  neutron  source  is 


anisotropic  because  of  its  shape  and  cladding.  The 
239pu-Be  source  is  a  cylinder  9.45  centimeters  long  and 
2.08  centimeters  in  diameter  enclosed  with  Tantalum  and 
stainless  steel.  Thus,  more  neutrons  will  exit  at  90°  to 
the  cylindrical  axis  than  axially  through  the  ends.  The 
anisotropy  factors  relative  to  isotropic  emission  for  the 
source,  M1170,  are  cited  as  1.071  at  90°  to  the  axis  and 
0.48  at  the  ends. 

A  four-mil  cadmium  sheet  was  placed  in  the  path  of  the 
neutrons  to  eliminate  the  thermal  neutron  contribution  to 
the  t ime-of-f 1 i ght  spectrum.  Cadmium  has  an  absorption 
cross  section  for  thermal  neutrons  of  about  19,910  barns, 
and  can  effectively  absorb  all  neutrons  below  0.5  eV . 

The  first  experiments  were  set-up  using  the  geometry 
in  figure  (3).  When  the  experiment  was  first  attempted, 
the  data  did  not  respond  as  predicted.  The  low  energy  end 
of  the  spectrum  was  very  high  and  the  upper  end  was  low. 
Investigation  into  this  revealed  that  the  second  detector 
was  receiving  over  a  factor  of  ten  times  the  expected 
number  of  counts.  To  reduce  this  number  and  bring  the 
counts  down  on  the  second  detector,  the  shadow  shield  was 
bu i It. 
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Material  for  the  shadow  shield  consisted  of  reactor 


grade  graphite  and  paraffin,  which  contains  a  large  amount 
of  hydrogen  atoms.  To  compute  the  thickness  of  graphite 
necessary  for  the  moderation  of  the  most  energetic 
neutron,  the  macroscopic  cross  section  for  the  atoms  was 
needed.  This  property  for  graphite  is  given  in  Table  1-4 
of  Nuclear  Reactor  Theory  by  LaMarsh  as  0.3851  cm  1  which 
equates  to  a  mean  free  path  for  the  neutron  to  be  2.59  cm. 
A  neutron  of  11  MeV  which  maintained  7 1 . 6^6  of  its  energy 
following  a  collision  with  a  carbon  atom,  would  require  a 
total  of  10  collisions  to  reduce  its  energy  down  to  less 
than  0.5  MeV  (7:168).  Therefore,  a  total  of  25.9 
centimeters  or  10.2  inches  of  graphite  were  needed  to 
moderate  the  neutrons  scattered  from  the  walls  and  floor. 
The  paraffin  was  added  to  ensure  the  moderation  of  the 
neutrons  prior  to  striking  the  detector.  The  detector  was 
covered  with  a  4-mil  sheet  of  cadmium  to  prevent  these 
thermal ized  neutrons  from  entering  the  detector.  The 
paraffin  and  carbon  were  not  placed  in  front  of  the 
sensitive  volume  of  the  detector  to  allow  the  scattered 
neutrons  from  the  first  detector  to  enter  the  sensitive 
volume  without  moderation.  With  the  shield  in  place,  the 
number  of  counts  on  the  second  detector  reduced  by  a 
factor  of  seven  over  that  of  an  unshielded  detector. 
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The  first  experiment  continued  for  104  hours  before 
the  time  spectrum  evolved  into  a  useful  piece  of  data. 

The  different  sections  of  the  spectrum  were  smoothed  by 
using  a  feature  on  the  dual  parameter  analyzer  known  as 
the  smooth  command.  This  feature  statistically  smooths 
i rregu 1 ar i t i es  in  spectra  without  altering  the  sizes  and 
shapes  of  the  spectral  peaks.  It  performs  a  5-point  data 
smoothing  of  the  spectrum  in  a  region  of  interest. 

Smoothed  data  (S)  are  calculated  from  spectral  data  (Y)  in 
channel  n  by  the  equation 


S(rc)  = 


WY(n)>J2(.Y(n-  1)*  Y(n+  1))-  3 (.Yin  2)*  Y(n  +  2)) 

" .  "  .  35  . . 


(24) 


Smoothed  data  allow  for  a  smoother  energy  spectrum  from 
the  unfolding  of  the  time  spectrum. 

In  the  second  measurement,  the  second  detector  was 
placed  65  centimeters  from  the  first  detector  and  the 
delay  time  in  the  anode  circuit  of  the  first  detector 
changed  to  100  nanoseconds.  The  calibration  and  timing 
were  repeated  using  the  procedure  described  in  Appendix  B. 
Experiment  2  produced  a  usable  spectrum  in  one  day; 
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however,  the  data  was  allowed  to  accumulate  over  the  same 
time  interval  as  experiment  1  to  check  the  number  of  data 
points  collected  for  the  same  time  interval. 

Experiments  3-6  were  conducted  in  a  low-scatter 
environment  with  no  shield  in  place.  All  of  these 
experiments  were  conducted  in  the  dome  of  building  470. 

The  geometry  of  the  experiment  is  shown  in  figure  4.  The 
scintillators  for  experiments  3  and  4  were  the  same  as 
those  used  in  experiments  1  and  2,  with  l-  1.5  meters  for 
experiment  3,  and  0.65  meters  for  experiment  4. 
Scintillator  selection  for  experiments  3  and  4  were  the 
same  as  those  used  in  the  first  two  experiments.  For  the 
last  two  experiments,  the  BICRON  422  scintillators  which 
measured  2"  in  diameter  and  1"  long  were  replaced  by 
BICRON  404  scintillators  which  measured  2"  in  diameter  and 
were  2"  long. 

The  distance  from  the  source  to  the  center  of  the 
sensitive  volume  of  the  first  detector  for  experiments  3 
and  4  was  34  centimeters.  Distances  between  centers  of 
the  sensitive  volumes  of  the  two  detectors  were  placed  at 
1.5  meters  for  experiment  3  and  5  and  0.65  meters  for 
experiments  4  and  6.  Several  days  lapsed  before  the 
spectrum  generated  from  each  experiment  could  be  used  to 
produce  an  energy  spectrum.  Once  the  time  spectrum  was 
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To  Crane 


Figure  5.  Experiment  3-6  Geometry 


generated,  it  was  once  again  smoothed  using  the  same 
procedure  previously  identified.  Experiment  4  was 
identical  to  experiment  2  with  the  exception  of  the 
location  of  the  experiment  itself. 
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Following  the  completion  of  experiment  4,  the 
scintillator  detectors  were  exchanged  for  the  2"  X  2" 
plastic  scintillators  with  the  titanium  oxide  coating. 

Experiments  5  and  6  were  the  same  as  experiments  3  and 
4,  with  the  exception  of  the  new  scintillators.  The  data 
for  experiment  5  did  not  produce  a  spectrum  which  could  be 
unfolded.  Analysis  determined  the  reason  for  this  is  the 
spread  in  the  angle  0.  With  the  first  detector  at  34 
centimeters  from  the  source,  this  spread  was  ±6°  which 
induced  an  error  in  the  measured  neutron  energy  of  =*20%. 
With  this  much  variation  in  the  energy  of  the  neutron,  the 
energy  spectrum  overlapped  itself,  rendering  it  useless. 

To  reduce  the  spread  in  scatter  angle,  the  first  detector 
was  moved  to  a  distance  of  1  meter  from  the  source.  This 
reduced  the  spread  of  the  angle  0  to  ±3°  and  the  error  for 
the  energy  reduced  to  less  than  ±10%. 

To  investigate  the  relationship  between  the  amount  of 
energy  deposited  in  the  recoil  protons  of  the  first 
detector  to  the  t ime-of - f 1 i ght  data  a  dual  parameter 
experiment  was  conducted.  Experiment  7  investigated  the 
dual  parameter  analysis  of  the  information  being 
generated.  The  proton  recoil  spectrum  from  the  first 
detector  was  put  into  the  y  axis  of  the  analyzer  and  the 
time-of-f 1 ight  data  into  the  x  axis.  Figure  (6)  is  a 
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picture  of  the  proton  recoil  spectrum.  The  time  of  flight 
data,  when  coincident  with  the  proton  recoil  data,  would 
generate  a  count  in  the  z  direction  of  the  analyzer.  The 
belief  was  that  the  t  ime-of-f  1  i  ght  data  was  prop:,  tional 
to  the  energy  deposited  in  the  first  detector.  Should  a 
low  energy  neutron  deposit  half  of  it’s  energy  into  the 
first  detector,  and  scattered  into  and  reacted  with  the 
second  detector,  then  the  small  pulse  height  generated 
should  be  in  coincidence  with  the  longer  t ime-of - f 1 i ght . 
The  expected  result  was  the  build  up  of  the  t ime-of -f 1 i ght 
data  along  a  diagonal  of  the  matrix.  This  experiment  was 
allowed  to  run  for  one  week.  The  results  are  investigated 
in  the  analysis  and  results  portion  of  the  report. 
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III.  Analysis  and  Results 


The  analysis  of  the  data  consists  of  several  phases. 
First  the  energy  spectra  of  the  different  experiments  are 
resolved  and  compared  to  the  other  spectra  from  other 
literature.  Next,  the  average  energy  of  the  spectrum  is 
calculated  by  determining  the  grand  mean  of  the  different 
spectra's  average  energy.  Finally,  an  error  analysis  is 
conducted  through  the  use  of  a  Monte  Carlo  simulation  to 
determine  the  energy  resolution  of  the  spectrum.  To 
ensure  the  linearity  of  the  T ime-to-Pu 1 se-He i ght 
converter,  the  calibration  was  conducted  prior  and 
following  each  experiment.  In  all  cases  the  value  of  the 
calibration  equation  did  not  change  from  the  start  of  the 
experiment  until  the  run  was  completed.  Therefore,  the 
assumption  of  the  linearity  of  the  circuitry  is  valid. 


Spectral  Data 

The  published  data  to  compare  the  energy  spectrum  from 
the  spectrometer  was  taken  out  of  an  article  published  by 
Arun  Kumar  and  P.  S.  Nagaran  in  1977  (8:176).  Figure  (7) 
is  a  picture  of  the  energy  spectrum  they  published.  From 
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this  normalized  spectrum,  the  peaks  at  3.1,  4.5,  6.5,  7,7, 
and  9.7  MeV  are  readily  distinguishable.  The  dip  at  6.1 
MeV  is  also  a  notable  feature  of  the  spectrum. 

In  the  high  scatter  environments  of  experiment  1,  the 
resolution  of  the  time  spectrum  took  a  week  to  develop. 

The  time  spectrum  was  slow  in  that  the  number  of  incident 
neutrons  on  the  second  detector  was  0.17  neutrons  per 
scatter  neutron  from  the  first  detector.  This  value  was 
determined  analytically  by  taking  the  solid  angle  of  the 
scatter  from  the  first  detector  to  the  second  detector  and 
determining  the  total  number  of  neutrons  incident  on  the 
second  detector  from  a  4n  scatter.  The  massive  scatter 
contribution  of  the  environment  where  the  experiment  was 
located  initially  generated  a  large  low-energy  spectrum 
which  did  not  match  the  spectrum  in  figure  (7).  The 
shadow  shield  helped  reduce  the  amount  of  scatter  and  the 
experiment  continued.  The  time  spectrum  from  the  analyzer 
is  shown  in  figure  (8).  The  times  for  the  different 
channel  numbers  are  obtained  by  use  of  the  calibration 
equation  for  the  experiment.  This  calibration  is  done 
using  the  procedure  described  in  Appendix  B.  Following 
the  calibration  the  analyzer  gives  a  slope  and  y  intercept 
for  the  calibration  equation.  A  typical  equation  for  the 
experimental  set-up  had  a  negative  slope  and  a  y  intercept 


49 


S(E)  MeV 


0.3 


Figure  7 


Published  Energy  Spectrum  of 
Source  by  Arun  Kumar  and  P.  S. 


a  5  Curie  239Pu-Be 
Nagara jan 
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close  to  the  value  of  the  delay  selected  for  the  signal 
from  the  first  detector.  The  time  for  each  channel  A t  , 
could  then  be  calculated  using  the  calibration  equation 


Channel  Number 


Figure  8.  Time  Spectrum  for  Experiment  1 
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A  TK  Solver  Plus  model  was  used  to  resolve  the  time 
spectrum  into  the  energy  spectrum.  This  model  consisted 
of  developing  a  normalized  spectrum  from  the  time  spectrum 
by  dividing  the  individual  components  of  the  spectrum 
(n(t))  by  the  integral  of  the  time  spectrum  from  0  to  °°. 
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With  the  normalized  components  of  the  time  spectrum 
S(A(),  the  values  of  the  energy  spectrum  could  be 
calculated  using  the  results  from  equation  (21).  Figure 
(9)  is  the  normalized  energy  spectrum  resolved  from 
experiments  1  and  2. 


As  depicted,  the  locations  of  the  peaks  and  valleys 
for  the  resolved  energy  spectrum  are  close  to  those  in 
figure  (8).  There  are  slight  differences  in  the  values 
for  the  fractional  amount  in  each  section  of  the  spectrum. 
These  differences  are  probably  caused  by  the  high  scatter 
environment  which  contaminated  some  portions  of  the  data 
resulting  in  this  poor  resolution. 

Experiments  3  and  4  resulted  in  a  better  energy 
resolution  spectrum,  due  to  the  elimination  of  most  of  the 
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Figure  9.  Normalized  Neutron  Energy  Spectrum  from 

Experiments  1  and  2 


scatter  contribution  from  the  concrete  walls.  Figures 
(10)  and  (11)  show  the  time  and  resolved  energy  spectra 
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produced  from  experiment  3.  The  peaks  for  the  different 
energy  resolution  in  the  time  spectrum  are  marked  for 
clarity. 


Figure  10.  Time  Spectrum  for  Experiment  3 
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Figure  11.  Neutron  Energy  Spectrum  for  Experiment  3 


Analysis  of  the  energy  spectrum  indicates  the  peaks 
and  valleys  match  the  published  literature  values.  The 
peak  at  3.2  MeV  is  of  the  appropriate  fractional  quantity; 


55 


however,  the  spectrum  still  continues  to  1 1  MeV  rathe^ 
than  10.5  MeV.  The  origin  of  this  error  is  still 
uncertain  but  believed  to  have  been  caused  by  the  arrival 
in  the  second  detector  of  scattered  neutrons  immediately 
after  the  excitation  of  the  first  detector  by  an  incident 
neutron.  Experiment  4  resulted  in  no  appreciable 
difference  in  the  resolved  energy  spectrum.  Following 
this  experiment,  the  scintillators  were  changed  and  the 
system  allowed  to  stabilize  for  twenty-four  hours. 

Experiment  5  indicated  that  the  number  of  interactions 
in  the  first  detector  increased  by  a  factor  of  2.5  over 
that  of  the  first  scintillators.  This  is  due  to  the 
doubling  of  the  hydrogen  atoms,  which  would  cause  the 
reaction  rate  of  the  incident  neutrons  to  double,  and  the 
doubling  of  the  path  length  which  the  neutrons  travel 
through  the  medium  which  increased  the  intrinsic 
efficiency  of  the  detectors.  Interactions  in  the  second 
detector  increased  by  a  factor  of  five.  After  two  days, 
the  time  spectrum  from  the  experiment  did  not  evolve  as 
the  ones  from  the  previous  experiments.  Investigation 
revealed  the  spread  in  the  scatter  angle  to  be  large. 

This  spread  was  from  0.7  radians  to  0.907  radians  which 
induced  a  large  error  in  the  resolution  of  the  energy 
spectra.  With  this  large  of  a  spread  in  0  ,  the  time 
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resolution  disappeared,  generating  a  spread  in  time  for  a 
5.25  MeV  recoil  neutron  to  vary  from  62.34  nanoseconds  to 
77.16  nanoseconds. 

To  correct  this,  the  first  detector  was  placed  1.0 
meters  away  from  the  center  of  the  source.  When  the 
experiment  was  repeated,  the  spread  in  the  scatter  angle 
and  time  reduced.  The  values  for  the  spread  in  scatter 
angle  went  from  the  values  previously  mentioned  to  0.75  to 
0.84  radians,  and  the  time  spread  reduced  to  65.0  to  71.0 
nanoseconds.  The  time  resolution  reappeared  and  the 
energy  spectrum  resolved  into  Figure  (12). 

To  compare  the  obtained  data  to  the  published 
spectrum,  the  spectrum  in  figure  6  was  over  laid  as  shown 
in  Figure  (13).  The  plot  indicates  that  the  spectrum 
between  3  and  6  MeV  did  not  correspond  with  that  in  the 
literature.  The  dip  at  6.1  MeV  occurs  at  5.9  MeV,  and  the 
fraction  of  neutrons  between  2  and  3  MeV  are  not  of  the 
appropriate  values.  Neutron  distribution  between  6  and 
10.5  MeV,  however,  are  virtually  the  same  as  the  one 
published.  Experiment  6  demonstrated  the  same  spectrum 
as  experiment  5  with  no  improvement  over  the  resolved 
energy  spectrum.  This  error  is  caused  by  the  energy 
resolution  from  the  spread  in  the  scatter  angle  0,  and  the 
error  induced  by  the  scatter  of  neutrons  in  the 
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Figure  12.  Experiment  5  Neutron  Energy  spectrum 


surroundings . 

To  conduct  the  dual  parameter  experiment,  several 
factors  had  to  be  implemented.  The  first  was  to  generate 
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and  attempt  to  calibrate  the  proton  recoil  spectrum  from 
the  first  detector  and  place  them  in  coincidence  with  the 
t ime-of -f 1 i ght  data.  This  spectrum  is  seen  in  Figure  (6). 
Calibration  of  the  proton  recoil  spectrum  was  not 
accomplished  due  to  the  lack  of  a  mono-energetic  neutron 
source  below  10  MeV .  Therefore,  the  experiment  was 
conducted  without  calibration. 

The-set  up  on  the  dual  parameter  analyzer  was  as 
follows:  the  t ime-of-f 1 ight  data  was  put  into  the  first 

channel  of  the  analyzer,  and  the  proton  recoil  spectra  was 
placed  on  the  second  channel.  A  coincidence  window  of  0.2 
micro  seconds  was  used  to  generate  the  data.  The  grid  was 
set  up  to  be  a  64  X  64  matrix.  If  a  t ime-of - f 1 i ght  pulse 
should  occur  during  the  coincidence  window  of  a  proton 
recoil,  then  the  analyzer  would  generate  a  count  at  that 
point.  To  ensure  that  the  timing  was  correct  for  the 
proton  recoil  spectra,  the  signal  from  the  linear  amp  of 
the  first  detector  was  delayed  the  same  amount  of  time  as 
the  anode  signal  from  the  first  detector.  Figure  (14)  is 
a  picture  of  the  analyzer  screen  following  a  week  of  data 
co 1 1 ect i on . 
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Figure  14.  Experiment  7  Dual  Parameter  Analyzer  Display 
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To  understand  the  plot,  the  following  information  is 
needed.  From  the  upper  right  corner  of  the  plot,  the 
t ime-of-f 1 i ght  information  is  along  the  axis  going  to  the 
left  and  the  proton  recoil  spectra  is  on  the  right  axis. 

As  shown  in  this  picture,  the  more  energetic  neutrons 
produce  the  longer  flight  time  and  generate  the  first 
peaks  from  the  left.  The  coincidence  window  of  200 
nanoseconds  resulted  in  the  generation  of  a  lot  of 
spurious  pulses  in  the  lower  ranges  under  this  peak,  thus 
reducing  the  readability  of  the  plot.  The  t ime-of-f 1 i ght 
data  can  be  determined,  however,  by  looking  at  the  highest 
proton  recoil  energy  of  each  t ime-of-f 1 i ght  peak.  If  this 
is  accomplished,  the  spread  of  the  t ime-of-f 1 i ght  data 
corresponds  to  that  previously  predicted.  The  neutrons 
deposit  half  of  their  energy,  (Cosz(^)  ),  into  the  hydrogen 
atom.  Therefore  the  lower  the  energy  of  the  neutron 
incident  on  the  first  detector,  the  longer  the 
time-of-f 1 ight  would  be.  Realizing  the  time  spectrum  is 
backwards  for  the  t ime-of-f 1 ight  data,  i.e.  channel  1  = 

200  nanoseconds,  channel  256  =  0  nanoseconds,  the  plot 
starts  to  come  together.  Additional  studies,  however; 
need  to  be  accomplished  to  reduce  the  spurious  pulses  from 
the  t ime-of-f 1 i ght  axis. 
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Average  Energy 


To  compute  the  average  energy  of  the  different 
spectra,  the  first  moment  of  the  spectra  was  divided  by 
the  zeroth  moment. 


/' 


ES(E)dE 


aug 


/; 


S(E)dE 


(26) 


This  integration  was  accomplished  numerically,  and  the 
grand  mean  of  the  average  values  of  the  different  energy 
spectra  was  determined  to  be  4.6±0.1MeK  .  This  value  is 
the  same  value  obtained  from  the  spectrum  in  figure  (6). 

To  calculate  the  average  efficiency  for  the  detector 
system  the  average  path  length  for  the  neutron  travel  in 
the  different  detectors  was  taken  as  2.54  centimeters  for 
experiments  1-4,  and  5.08  centimeters  for  experiments  5 
and  6.  Using  the  value  for  the  average  energy  previously 
obtained,  and  equations  (22)  and  (23),  the  average 
efficiency  for  the  detectors  was  calculated  to  be  0.1  ±0.02 
for  experiments  1-4  and  0.2±0.02  for  experiments  5  and  6. 
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Error  Analysis 


The  error  analysis  for  the  experimentation  undertaken 
proved  to  be  very  difficult.  From  a  single  deviation  in 
the  amount  of  time,  length  of  travel,  and  scatter  angle, 
the  errors  needed  to  be  propagated  from  equation  (9) 
through  equation  (21).  Even  if  Poisson  statistics  are 
assumed,  the  propagation  of  the  various  errors  was 
insurmountable.  The  energy  resolution  is  the  main 
concern  and  to  determine  the  spread  in  this  value  the 
scatter  angle  spread  must  be  known. 

Defining  the  scatter  angle,  6  ,  consumed  a  large 
portion  of  the  analysis  phase  of  this  project.  For  the 
first  two  experiments,  the  spread  of  e  was  0.1  radians  and 
0.14  radians,  respectively.  These  values  were  obtained 
through  the  error  analysis  performed  by  using  a  Monte 
Carlo  simulation  routine. 

The  Monte  Carlo  simulation  routine  established  for  the 
error  analysis  was  constructed  through  the  use  of  a  TK 
Solver  model.  A  copy  of  this  model  is  listed  in  Appendix 
C.  This  model  took  1000  random  points  for  each  of  the 
following:  (1)  the  place  of  the  neutron  generation  in  the 

source;  (2)  place  for  the  scatter  in  the  first  detector; 
and  (3)  place  for  the  interaction  in  the  second  detector. 
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Then  for  the  same  time  of  flight.  At  ,  the  energy  of  the 
neutron  incident  on  the  first  detector  was  calculated. 
Accomplishments  of  this  was  undertaken  by  generating 
non-nominal  distances  of  flight  and  scatter  angle,  e  . 

With  these  values,  the  energy  of  the  scattered  neutron 
could  be  calculated  based  on  the  distance  traveled  and  the 
time  it  took  to  traverse  this  distance  (constant).  The 
scatter  angle  determined  for  this  event,  and  the  energy  of 
the  incident  neutron  calculated  through  the  use  of 
equation  (22).  The  details  of  the  sampling  method 
employed  in  the  error  analysis  model  follows. 

To  sample  the  probability  of  the  location  of  the 
different  events  the  use  of  a  probability  distribution 
function  was  needed.  In  order  to  generate  this  function, 
the  probability  of  an  interaction  was  taken  to  be  that  of 
a  differential  area  of  a  cylinder, 

P(r  ,0 ,  z^drdQdz  =  rdrdQdz  (27) 

Which  indicates  that  the  probability,  P,  must  be 
proportional  to  the  location  of  the  event  in  the  radius, 
r . 

To  sample  the  event,  inversion  of  the  probability 
distribution  function  was  necessary.  This  was 
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accomplished  by  taking  the  ratio  of  the  respective 
circular  areas  for  the  location  of  the  event  at  radius,  r, 
to  that  of  the  circular  area  of  the  entire  radius,  R,  and 
setting  this  value  equal  to  a  random  number  £. 


or  simply. 


(28) 


(29) 


Therefore,  by  use  of  cylindrical  geometry,  the  value 
of  the  sampled  radius  of  the  event,  r,  the  angle  e  ,  and 
height,  z,  are  given  by: 


r-Ryfc 

z  =  h%  (30) 

0  =  2n% 

where  R  is  the  radius  of  the  different  components,  and  h 
is  the  height  of  the  different  components. 

The  value  of  £  for  the  z  component  is  computed  a  little 
differently  than  just  from  the  generation  of  a  random 
number.  Because  the  flux  of  neutrons  will  decay 
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exponentially  through  the  sensitive  volume  of  the 
detector,  the  random  number  must  be  biased  toward  the 
front  of  the  detector.  To  accomplish  this  task,  the 
cumulative  distribution  function 

C.  D.  F .  =  ( 1  -exp(-ad))  (31) 

is  used  to  determine  the  value  of  interaction  within  the 
volume  d,  where  d  is  generated  by  taking  the  negative  of  a 
natural  log  of  a  random  number  divided  by  the  value  of  the 
macroscopic  cross  section  a.  Once  this  value  is  known 
the  value  of  £  is  determined. 

With  the  different  points  established  for  the 
different  events,  the  analysis  of  the  energy  resolution 
could  begin.  To  accomplish  this,  the  initial  data  for  the 
six  experiments  was  placed  in  the  Monte  Carlo  routine. 

The  times  selected  for  analysis  corresponded  to  the 
fastest  neutron,  10.5  MeV,  the  slowest  neutron,  0.5  MeV 
and  the  average  neutron,  4.6  MeV.  The  data  for  the 
different  energies  for  these  times  were  then  sorted  using 
a  quick  basic  routine  and  the  number  per  group  plotted. 
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Figure  15.  Error  Analysis  for  High  Energy  Spectrum 
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Figure  16.  Er>~or  Analysis  for  Mid-Range  Neutron  Fnergy 


Analysis  of  the  data  in  figures  (15-17)  indicates  that 
the  average  energy  of  the  high  energy  spectrum  A t  ,  is 
given  as  1080  <  0.34  MeV  ,  for  the  mid  range  as  3.4  ±  0.04 

and  for  the  low  end  1.0/  >  o.oub  .  A  similar  analysis  was 
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Figure  17.  Error  A-a'ysis  r  . ow  Energy  Neutron 


accomplished  for  the  value  of  0  ,  and  length  I  .  It  was 

assumed  that  0  =  45°  for  the  analysis.  The  values 
obtained  for  the  different  components  were  0  -  45. 2C  ±-  0.2 

and  /  1.51  *  0.01  meters. 
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IV.  Discussion.  Conclusion s  and  Recommendat i ons  for  Futu r e 

Studies 


Energy  spectrum  from  a  239Pu-Be  source  was  determined  by 
using  a  neutron  t ime-of-f 1 i ght  spectrometer.  This 
spectrometer  measures  the  time  of  flight  of  the  neutron 
between  two  points  and  generates  a  corresponding  time 
spectrum.  The  time  spectrum  was  unfolded  into  the  energy  of 
the  incident  neutron  through  the  use  of  a  TK  Solver  model. 
Energy  spectrum  which  were  generated  from  the  time  spectrum 
were  compared  to  those  in  the  literature.  Several  experiments 
were  undertaken  to  investigate  the  effect  that  the  distance 
between  two  detectors,  l  ,  had  on  the  energy  resolution  along 
with  the  effect  of  two  different  scintillators.  The  two 
values  of  l  investigated  were  0.65  meters  and  1.5  meters. 
Scintillators  of  2"  in  height  were  replaced  for  ones  with  1" 
height  and  the  different  energy  spectra  considered. 
Additionally,  the  effect  of  the  high-scatter  environment 
contributed  to  the  building  of  a  shadow  shield  in  order  to 
prevent  the  scattered  neutrons  from  the  walls  and  floors  to 
corrupt  the  data.  Elimination  of  this  nuisance  was 
accomplished  by  moving  to  a  ' . w  scatter  environment. 

Uncertainty  in  the  time  spectrum  can  be  from  a  number  of 
sources.  Definition  of  the  scatter  angle,  0  ,  along  with  the 


determination  of  the  true  path  length  of  flight  are  the 
largest  contributors.  Multiple  scatters  of  the  neutron  with 
hydrogen  or  even  the  carbon  atoms  can  also  contribute  to  the 
uncertainty  in  the  unfolded  energy  spectrum. 

The  last  experiment  considered  the  coincidence  of  the 
proton  recoil  energy  spectrum  to  that  of  the  t ime-of - f 1 i ght . 
Even  though  the  results  are  hazy,  the  desired  effect  can  be 
seen.  Additional  study  should  be  undertaken  to  resolve  the 
long  coincidence  window  timing  problem. 

Although  several  types  of  spectrometers  have  been 
developed  to  determine  the  energy  spectra  of  different  neutron 
sources,  most  have  limitations  as  to  the  range  of  energy 
coverage.  The  t ime-of -f 1 i ght  spectrometer  has  emerged  as  the 
preferential  choice  for  determining  the  energy  spectrum  from 
any  neutron  source.  From  the  time  spectrum,  which  is 
generated  by  the  t ime-of - f 1 i ght  spectrometer ,  the  energy 
spectrum  is  easily  achieved  through  the  use  of  manipulation  of 
the  kinematics  equations  and  distribution  functions.  The 
error  analysis  of  the  energy  resolution  which  is  not  readily 
obtainable  from  the  propagation  of  errors  through  the 
equations,  can  be  accomplished  by  use  of  a  Monte  Carlo 
simulation. 

Neutron  spectrum  obtained  from  a  5  curie  Pu-Be  source 
with  the  t ime-of -f 1 i ght  neutron  spectrometer  are  approximately 
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as  described  in  the  literature.  They  take  a  long  time  to 
obtain  with  "weak"  sources  when  just  two  detectors  are  used. 

A  low-scatter  environment  is  essential  although  some  results 
are  possible  in  a  high-scatter  environment  if  appropriate 
shielding  is  used.  Refinements  are  required  to  improve  the 
resolution  of  the  time  spectrum,  such  as  defining  the 
dependence  of  0  on  the  different  geometries  available. 

Another  avenue  to  investigate  would  be  the  use  of  a  small 
scintillator  on  the  first  detector,  and  a  large  scintillator 
on  the  second  detector.  This  would  allow  the  placement  of  the 
first  detector  closer  to  the  source,  increasing  the  number  of 
neutrons  incident  on  it,  without  loosing  the  definition  of  the 
scatter  angle.  The  larger  scintillator  on  the  second  detector 
would  improve  the  efficiency  of  the  second  detector  and  reduce 
the  amount  of  time  necessary  to  generate  the  time  spectrum. 
Again  consideration  must  be  given  to  the  scatter  angle  spread. 

The  average  energy  for  the  neutron  source  selected  is 
46*0.1  Me V  which  is  the  same  value  from  the  spectrum  given 
by  Arun  Kumar  and  P.  S.  Nagarajan.  With  this  value  for  the 
average  energy,  the  intrinsic  efficiency  of  the  detectors  is 
determined  to  be  around  0.1  -  0.2.  This  efficiency  increases 
with  an  increase  in  the  volume  of  the  sensitive  volume  and  an 
increase  in  the  average  path  length  of  travel. 
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Further  studies  should  be  undertaken  to  improve  the 


resolution  of  the  time  spectrum,  and  thus  the  energy  spectrum 
of  the  source.  One  modification  could  be  the  placement  of  the 
second  detector  at  a  distance  of  3.0  meters  from  the  first 
detector  with  the  first  detector  at  a  distance  of  2.0  meters 
from  the  source.  This  would  allow  for  the  spread  in  the 
scatter  angle,  6  ,  to  be  minimized.  A  second  study  should  be 
undertaken  to  refine  the  dual  parameter  aspect.  By  placing  a 
greater  distance  between  the  detectors,  the  resolution  of  the 
coincidence  events  could  be  better  understood.  The  minimum 
coincidence  window  for  the  dual  parameter  analyzer  is  0.1 
microseconds  and  the  average  time  of  flight  for  the  neutron 
between  the  detector  should  be  at  least  an  order  of  magnitude 
greater  than  this.  The  use  of  several  detectors  in  an  array 
a  fixed  distance  from  the  first  detector  which  would  generate 
the  same  scatter  angle  would  reduce  the  amount  of  time 
necessary  to  generate  the  time  spectrum. 
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Appendix  A.  _ Equipment  list 

The  following  is  a  list  of  equipment  used  to  establish  the 
spectrometer  circuit  for  the  t ime-of -f 1 i ght  system. 


ITEM 

MANUFACTURER 

MODEL  NUMBER 

ID  NUMBER 

Photomultipl ier 

RCA 

8575 

P52180 

T  ube 

T30499 

Photomu 1 1 ip  1 ier 

Ortec 

265 

1589 

Base 

1592 

P  1  ast i c 

B  i  cron 

422  (Exp  1-4) 

none 

Scinti 1 lator 

404  (Exp  5-7) 

none 

Preamps 

Ortec 

113 

6511 

6571 

Linear  Amps 

Ortec 

572 

0978 

4745 

Quad  CFD 

T  enne 1 ec 

TC  454 

509 

i  Gate  and  Delay 

Ortec 

4  1  6  A 

3640 

Gener ators 

3639 

3679 

1 


ITEM 

MANUFACTURER 

MODEL  NUMBER 

ID  NUMBER 

Timing  SCA 

Ortec 

551 

2755 

2768 

T imer/Counter 

Ortec 

996 

116 

125 

Pulse  Shape 

Analyzer 

Ortec 

458 

567 

Delay  Amp 

Ortec 

463 

832 

Time  to  Amplitude 

Converter 

Tennelec 

TC-873 

244 

Fast  Coincidence 

Circuit 

Ortec 

414A 

810 

Time  Cali brator 

Ortec 

462 

390 

High  Voltage 

Power  Supply 

Ortec 

556 

1993 

2 


ITEM  MANUFACTURER  MODEL  NUMBER  ID  NUMBER 


Osci 1 1 oscope 

T  ektron i cs 

425 

3465 

NIM  Bins 

Ortec 

400 1C 

3583 

3520 

Dual  Parameter 

Nuclear  Data 

ND  76 

963085 

Analyzer 

Neutron  Source 

Monsanto 

Pu-Be 

Ml  170 

3 


The  following  procedure  was  used  to  time  the  circuitry  and 
the  calibrate  the  time  to  amplitude  converter  for  use  in  the 
spectrometer . 

1.  Place  both  detectors  equal  distance  from  neutron  source. 

2.  Feed  inputs  to  the  fast  coincidence  circuit  from  detector 

2  to  the  dual  channel  oscilloscope.  Select  triggering  mode  to 
channel  one  or  channel  2. 

3.  Adjust  gate  and  delay  circuits  until  pulses  are  in 
coincidence  on  scope. 

4.  Re-attach  cables  to  inputs  of  fast  coincidence  circuit; 
set  resolving  time  to  100  nanoseconds. 

5.  Check  output  of  fast  coincidence  circuit;  should  be 
positive  logic  pulse.  If  not  repeat  steps  2-4. 

6.  Remove  all  delays  from  detector  1  anode  circuit. 

7.  Place  start  and  stop  signals  from  the  time  to  amplitude 
converter  to  the  dual  channel  oscope. 

8.  Adjust  gate  and  delay  circuit  of  the  anode  signal  from 
detector  1  only  to  bring  signals  into  coincidence;  trigger 
oscope  on  start  signal  from  detector  2. 

9.  Insert  fixed  delay  into  anode  signal  from  detector  1. 

10.  Re-attach  start  and  stop  signals  to  time  to  amplitude 
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converter . 


11.  Select  range  scale  on  time  to  amplitude  converter  to 
match  fixed  delay  set  in  step  9. 

12.  Place  multiplication  plot  to  1.25. 

13.  Check  to  see  that  you  have  valid  conv  lite;  tac  out  lite; 
and  SCA  lite  on . 

14.  Remove  the  start  and  stop  signals  from  the  time  to 
amplitude  converter . 

15.  Place  the  start  and  stop  signals  from  the  time  calibrator 
into  the  start  and  stop  in  connections  of  the  time  to 
amplitude  converter . 

16.  Select  range  on  the  time  to  amplitude  converters  to  a 
value  just  under  the  fixed  delay  selected  for  the  anode 
circuit  from  detector  1. 

17.  Select  period  setting  to  20  nanoseconds. 

18.  Acquire  data  from  time  to  amplitude  converter  on  dual 
parameter  analyzer. 

19.  Using  RM  and  LM  commands  place  right  and  left  markers  on 
the  peaks  generated  by  the  time  to  amplitude  converter. 

20.  Determine  the  value  of  the  right  peak  and  left  peak  by 
subtracting  the  value  of  the  range  determined  in  step  16  from 
the  fixed  delay  in  the  anode  circuit  of  detector  1;  this  is 
the  value  of  the  right  marker  energy. 

21.  The  other  peaks  will  be  every  20  nanoseconds  apart, 
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increasing  from  right  to  left,  instead  of  left  to  right. 

22.  Use  the  calibration  commands  LE  and  Rl  and  place  the 
value  of  the  peaks  for  the  left  and  right  markers  into  the 
dual  parameter  analyzer. 

23.  Remove  the  start  and  stop  signals  from  the  time 
calibrator  and  reconnect  start  and  stop  signals  from 
spectrometer  circuit. 

24.  Reestablish  geometry  for  experiment  and  commence 
experiment . 
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Appendix  C. _ TK  Solver  Models  Used  for  Analysis 


The  following  model  is  the  one  used  to  resolve  the  time 
spectrum  into  the  energy  spectrum.  To  accomplish  this  the 
model  first  generates  a  normalized  time  spectrum  and  then  uses 
equation  17  to  generate  the  energy  spectrum.  This  spectrum  is 
then  analyzed  for  the  average  energy  by  taking  the  value  of 
the  first  moment  of  the  spectrum  over  the  zeroth  moment  of  the 
spectrum.  The  function  for  integration  is  a  list  function 
with  cubic  interpolation  between  the  data  points. 


Variable  Sheet 


St  1 nput 

Name 

Output 

Unit 

Comment 

E 

Jou 1 es 

recoil  neutron  energy 

Er 

MeV 

recoil  neutron  energy 

1 . 627E-27 

m 

Kg 

mass  of  neutron 

V 

m/sec 

velocity  of  neutron 

1.5 

1 

meters 

distance  between  det 

1  and  det  2 

del  tat 

sec 

time  of  travel 

0 

Ch 

channel  number 

C 


1 


st 

0  nt 

0  Se 

0  En 

45  theta 

Num 

0  1  ow 

11  high 


1000  div 

Ru  1  e - 

*  Er  =  E/1 .602e-13 

*  E  =  . 5*m*v" 2 

*  v  =  1/deltat 


sec''  - 1 

MeV~-1 

MeV 

degrees 

.12722 


normal i zed  time 
spectrum 
time  spectrum 
Energy  spectrum 
Incident  Neutron 
Energy 

Scatter  angle 
Value  of  first  moment 
smallest  value  of 
energy 

highest  value  of 
energy 

Zeroth  moment 
Average  Energy  of 
spectrum 

Number  of  bins  for 
numerical  integration 


Den  .035327 

Eavg  3.60124 


Converts  from  Joules  to  MeV" 
"Kinematics  Equation  for 
Kinetic  Energy" 

"Definition  of  velocity" 
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*  del  tat  =  ( - . 8*Ch+ 196.8)*e-9 

*  st  =  nt/sum( ' nt ) 

*  Se  =  st*(deltat~3)/(m*l ^2) 

*  Er  =  (cosd(theta) )  ,'2*En 

*  Num  Fmom( low, high) 

*  Den  Trap( low, high) 

*  Eavg  =  Num/Den 


"Calibration  equation  for 
exper iment 

"Determines  s(t)  from  n(t) 
"Determines  s(e)  from  s(t)" 
"Determines  neutron  energy  from 
recoil  energy" 

"Determines  first  moment  for 
spectra  using  trapezoid  rule" 
"Determines  zeroth  moment  for 
spectra  using  trapezoid  rule" 
"Determines  avg  energy" 


The  second  model  used  is  the  error  analysis  model.  To 
determine  the  error  in  the  energy  spectrum,  the  random  number, 
|  ,  is  used  along  with  equation  26.  The  zero  plane  for  the  z 
axis  is  the  bottom  of  the  neutron  source,  whereas  the  zero 
plane  for  the  y  axis  is  the  bottom  of  the  first  detector's 
sensitive  volume  and  the  zero  plane  for  the  x  axis  is  the 
center  line  of  the  first  detector's  sensitive  volume. 

Using  this  for  the  origins  we  see  that  the  values  for  the 
location  of  the  neutron  generation  in  the  source  (x,  y,  z)  are 
given  by  (rcos  0,  rsin  0,  z)  .  The  values  for  the  neutron 
reaction  (x,  y,  z)  in  the  first  detector  are  given  by  (rcos 
0,  z,  rsin  0).  The  last  detector  has  initial  coordinates  based 
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on  the  geometry  of  the  experiment  and  then  the  values  of  the 
neutron  interaction  (x,  y,  z)  are  given  by  the  initial  value 
of  the  bottom  center  of  the  sensitive  volume  minus  (zcos  - 
rcos  e  sin  rsin  0).  Values  for  the  time  of  travel,  energy 
of  neutron,  and  spread  in  0  can  be  determined.  It  is 
important  to  note  that  you  must  input  either  the  energy  of  the 
neutron  in  order  to  generate  a  time  difference  or  a  fixed  time 
to  generate  the  energy  spread. 

=  VARIABLE  SHEET  ■  . . , For  Academic  Use  Only 

St  Input -  Name -  Output -  Unit -  Corrment - 

INITIAL  PARAMETERS  FOR  PROBLEM 

************************************* 


150 

LDET1TO 

cm 

Distance  from  detl  to  det2 

100 

LSORTOD 

cm 

Distance  from  source  to  detl 

1.0414 

Rsource 

cm 

Radius  of  source 

9.488 

Zsource 

cm 

height  of  source 

2.54 

Rdetl 

cm 

radius  of  detl 

2.54 

Zdetl 

cm 

length  of  det  1 

2.54 

Rdet2 

cm 

radius  for  det  2 

2.54 

Zdet2 

cm 

length  of  det  2 

46.904442 

del tat 

nanosec 

Time  for  travel  between  detl 

det2 

1.675E-27 

M 

kg 

mass  of  neutron 
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GENERATE  RANDOM  VALUES 


************************************* 


L  .2 

random 

random  number 

rsource 

cm 

random  r  for  source 

zsource 

cm 

random  z  for  source 

thetaso 

radians 

random  single  for  source 

rdetl 

cm 

random  radius  for  detl 

zdetl 

cm 

random  z  for  det  1 

thetade 

radians 

random  angle  for  det  1 

rdet2 

cm 

random  radius  for  det 2 

zdet2 

cm 

random  z  for  det  2 

thetade 

radians 

random  angle  for  det  2 

OONVERT  TO  RECTANGULAR  COORDINATES 


************************************* 


xsource 

cm 

X 

component  of 

source 

ysource 

cm 

y 

component  of 

source 

ZSOURCE 

cm 

z 

component  of 

source 
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xdetl 

cm 

x  component  of  detl 

ydetl 

cm 

y  component  of  detl 

ZDET1 

cm 

z  component  of  detl 

xdet2 

cm 

x  component  of  det2 

ydet2 

cm 

y  component  of  det2 

ZDET2 

cm 

z  component  of  det2 

ANALYS 1 S 

************************************** 

CALCULATED  VALUES  FOR  ERROR 

lendetl 

cm 

Length  detl  to  det2 

L  .76492326  theta 

radians 

for  neutron  travel 

Scatter  angle 

El 

MeV 

recoil  energy  of  neutron 

L  10.168283  En 

MeV 

Incident  neutron  Energy 

************************************* 

MINIMUM  AND  MAXIMUM  VALUES 

FOR  ANALYSIS 

minimum 

MeV 

Minimun  value  of  En 

maximum 

MeV 

Maximum  value  of  En 

c 


6 


Minthet 

Maxthet 


radians  Minimum  theta(scatter  angle) 
radians  Maximum  theta (scatter  angle) 


mini eng 
maxi  eng 
mintime 
maxtime 


cm 

cm 

nanosec 

nanosec 


minimum  length  travel 
maximum  length  travel 
minimum  time  for  travel 
maximum  time  for  travel 


.78875847  thetaav  radians  average  theta 

L  thetasp  radians "2  intermediate  term  needed  for 

spread 

thetava  .00042376  radians''?  variation  of  theta 


L 


10.84695 


nbar 

MeV 

average  of 

data 

ssqddat 

MeV~2 

data  for  s' 

2 

ssqd 

MeV  "2 

s"2 

intermediate  values  necessary  for 
generation  of  data 


************************************* 


v 

thetaso 

thetade 


m/sec  velocity  of  recoil  neutron 
angle  source  to  det  1 
angle  detl  to  det  2 
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seed2 

random  #  generated 

seed3 

random  #  generated 

seed4 

random  #  generated 

seed5 

random  #  generated 

seed6 

random  #  generated 

seed7 

random  #  generated 

seed8 

random  #  generated 

seed9 

random  #  generated 

seed 10 

random  #  generated 

total E 

MeV 

10.168283  Eavgl 

barns 

Average  Energy  of  spectrim 

.91257036  signasl 

MeV 

microscopic  cross  section  for 

Hydrogen 

.04754492  Signasl 

cm"- 1 

Macroscopic  cross  section 

21.032743  Lambda  1 

cm 

atomic  density  for  hydrogen 

1.521548  sigmas 

barns 

microscopic  cross  section  for 

Hydrogen 

5.2922779  Eavg 

MeV 

Average  Energy  of  spectrum 

.07927265  Signas 

cm  -1 

Macroscopic  cross  section 

5.21E22  Nh 

atoms/cm" 

atomic  density  for  hydrogen 

12.614691  Lambda 

cm 

mean  free  path  for  det2 
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=  RULE  SHEET 


For  Academic  Use  Only 


S  Rule - 

seed2=SEED1 (random) 
r sour ce=Rsour ce* ( X i ( seed2 ) ) ~ . 5 
seed3=SEED1 (seed2) 
zsource=Zsource  *  Xi(seed3) 
seed4=SEED1 ( seed3) 
thetasource=2*pi ( )*seed4 

seed5=SEED1 (seed4) 
r det 1 =Rdet 1 * ( X i ( seed5 )  )  ~ .  5 
seed6=SEED1 (seed5) 
zdet1=Zdet1*zdet(seed6,Signas1 ) 
seed7 =SEED 1 ( seed6 ) 
thetadet1=2*pi ( )*seed7 

seed8=SEED1 (seed7 ) 
r det 2=Rdet 2* ( X i ( seed8 ) ) A . 5 
seed9=SEED1 (seed8) 
zdet2=Zdet2*zdet ( seed9 , S i  gras ) 
seed10=SEED1 (seed9) 
thetadet2=2*pi ()*seed10 

xsour ce=r sour ce*cos ( thetasource ) 
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ysource=LSORTODET 1 +rsource*s i n( thetasource ) 

ZSOURC£=z source 

xdet 1 =rdet 1 *cos ( thetadet 1 ) 

ydet1=zdet1 

ZDET1=. 06908  +  rdet 1*sin( thetadet 1 ) 

xdet2=(LDET1TODET2+Zdet2)*cos(-pi( )/4)  -  zdet2*cos(-pi ( )/4) 
ydet2=(LDET1TODET2+Zdet2)*sin(-pi ( )/4)  +rdet2*cos(thetadet2)*sin(pi ( )/4) 
ZDET2=. 06908  +  rdet2*sin(thetadet2) 


thetadet 1 todet 2=ABS ( atan2 ( ydet 1 +ydet 2 , xdet 1 +xdet2 ) ) 
thetasoutodet 1 = ABS ( atan2 ( ydet 1 +ysource , xdet 1 +xsource ) ) 
theta=p i ( ) - ( thetasoutodet 1 +thetadet 1 todet2 ) 

1 endet 1 todet2= ( ( xdet 1 +xdet 2 ) “ 2+ ( ydet 1 +ydet 2 ) * 2+ ( ZDET 1 +ZDET2 ) ~  2 ) ~ . 5 

v= 1 endet 1 todet 2/de 1 tat 

E1=.5  *  M*v^2 

En=E1/(cos(theta)',2) 

minimLir£=MIN(  'En) 

max  irmur£=MAX(  ’  En) 

Mintheta=MIN( ' theta) 

Maxtheta=MAX( 'theta) 

mini ength=M I N ( ' 1 endet 1 todet  2 ) 
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max1ength=Max( ' lendet1todet2) 
mintime  =  Min('deltat) 
maxtime  =  Max('deltat) 
tota1E=Sun( ’En) 
nbar =tota 1 E/ 1 000 
ssqddata= (En-nbsr ) " 2 
ssqd=Sum( ’ssqddata)/(999) 

*  Eavg 1 =En/ 1 . 602e- 1 3 

*  sigras1=4.83/((Eavg1)*.5)-.578 

*  Sigras1=Nh*signas1*1e-24 

*  Lambda1=  1/Sigrasl 

*  Eavg=E 1 / 1 . 602e- 1 3 

*  sigras=4.83/((Eavg)~ .5)- .578 

*  Sigmas=Nh*sigras*1e-24 

*  Lambda^  1 /Si gras 

*  thetaavg=Sum( ’theta)/ 1000 

*  thetaspread=(theta-thetaavg),'2 

*  thetavar=Sun( * thetaspread)/999 


—  PROCEDURE  FUNCTION:  Xi  ==» - =  ■  For  Academic  Use  Only 

Corrment:  Generates  random  numbers 


-  1 1 
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Parameter  Var i ab 1 es :  random 


Input  Variables:  seedl 

Output  Variables:  xi 

S  Statement - 

MaxLong:=  2147483647 
ModFactor : =7 . 826369259425646e-6 
if  seed 1=0  then  seed1=random 
y=seed1*MaxLong- . 5 
x=ModF actor  *  ROUND(y) 
seed2=x- INT(x) 
if  seed2  =  1  then  xi=random 
if  seed2<1  then  xi=seed2 
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1E24 
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Stop  on  List  Error:  No 


C 


12 


Use  Automatic  Iteration: 


Yes 


Comparison  Tolerance: 

.000001 

Typical  Value: 

1 

Maximum  Iteration  Count: 

20 

Global  Numeric  Format: 

Append  Variable  Names: 

Yes 

Use  Page  Breaks: 

No 

Number  Pages: 

Yes 

Form  Length: 

66 

Printed  Page  Length: 

60 

Printed  Page  Width: 

80 

Left  Margin: 

0 

Printer  Device  or  Filename: 

Printer  Setup  String: 

e :  \msdat\error 

Plot  Output  Filename: 

TKPLOT . OUT 

Use  Color: 

Yes 

Slow  Redisplay: 

Yes 

Solid  Line  Headings: 

Yes 

Bottom  Prompt  Line: 

Yes 

C  -  13 


vita 


Lieutenant  Commander  Wood  was  born  on 
Jackson,  Mississippi.  Shortly  after  graduation  from  Central 
High  School  in  1971,  he  joined  the  Navy  and  proceeded  to  serve 
in  several  combat  roles  in  Southeast  Asia.  In  1975  he  applied 
for  and  was  accepted  into  the  Navy  Enlisted  Scientific 
Education  Program  (NESEP),  where  he  earned  his  Bachelor  of 
Electrical  Engineering  degree  from  Auburn  University.  Upon 
graduation  he  was  commissioned  an  Ensign  in  the  United  States 
Navy  and  served  as  the  Assistant  Weapons  Officer  on  board  the 
U.  S.  S.  Lafayette  (SSBN  -  616  Blue)  and  as  the  Weapons 
Officer  on  board  the  U.  S.  S.  Von  Steuben  (SSBN  -  632  Blue) 
until  1985.  Following  his  tours  at  sea  he  changed  designators 
to  become  a  Engine**  : ng  Duty  Officer,  where  he  served  at  the 
Polaris  Missile  Far  •  1  :ty  Charleston,  South  Carolina.  While 
serving  in  this  billet,  he  was  responsible  for  development  of 
the  Limited  Life  Component  Exchange  (LLCE)  program  for  the 
Mark  IV  warhead.  Upon  activation  of  the  program,  he  entered 
the  School  of  Engineering,  Air  Force  Institute  of  Technology, 
in  August  1989. 


Permanent  Address: 


75 


References 


1.  Allen,  William  D.  Neutron  Detection.  New  York: 
Philosophical  Library  tnc,  1960. 

^  2.  A jzenberg-Se 1 ove ,  Fay.  Nuclear  Spectroscopy.  New  York: 

Academic  Press  Inc,  1960. 

3.  Kirk,  F.  W.,  "Neutron  T ime-Of -F 1 i ght  Spectrometers , " 

^  Nuclear  Instruments  and  Methods,  762:  539-563  (1979). 

4.  Knoll,  Glenn  F.  Radiation  Detection  and  Measurement.  New 

York:  John  Wiley  &  Sons,  1989. 

^  5.  Benedict,  Manson  And  Pigford,  Thomas  H.  And  Levi,  Hans  W. 

Nuclear  Chemical  Engineering.  New  York:  McGraw-Hill  Book 
Company ,  1981. 

*  6.  Cross,  William  G.  And  Ing,  Harry.  "Neutron  Spectroscopy," 

The  Dosimetry  of  Ionizing  Radiation,  Volume  2,  edited  by 
Kenneth  R  Kase.  Orlando:  Academic  Press  Inc,  1987. 

^  7.  LaMarsh,  John  R.  Introduction  to  Nuclear  Reactor  Theory. 

Reading  Massachusetts:  Add i son-Wes 1 ey  Publishing  Company, 
1972. 

*  8.  Kumar,  Arun  and  Nagarajan,  P.  "Neutron  Spectra  of  239Pu-Be 

Neutron  Sources,"  Nuclear  Instruments  and  Methods,  140: 
175-179  (1977). 

i 

> 


76 


;  REPORT  3 0 C U i'/l E MTATI G M 

PAGE 

1  . .  .  ■■>„*  V.’  • "«  -  -  -  .  .  .  .  ■  ■  .  „  •  -  “  -  -  3 

11  AGc.NC'  -ii  3 >i L IE.  REPORT  2A7E 

i  March  ] 99 ] 

j  i.  REPORT  7V  =  E  AND  DATES  COHERED 

Master’s  Thesis 

A.  TITLE  AND  EU3TI7U 

’  5.  EUNCJNG  Vu.V3  =  3S 

t 

•  Construction,  Calibration  and  Use  of  a 
Neutron  Time-of-Flight  Spectrometer 

i 

i 

i 

:  5.  AUTrlCRlS) 

i 

i 

!  Charles  L.  Wood  LCDR,  USN 

i 

! 

7-  ’E.RFCRMING  ORGANIZATION  NA.VIEiS)  AND  ACCRESSiESi 

Air  Force  Institute  Of  Technology 

WPAFB  OH  45433-6583 

•3.  ’E  REGRMISJu  ORGANIZATION 
•  REPCR'  IL.M3ER 

’  AFIT/GNE/ENP/9 ]M- ] 2 

:•  eronscring  monitoring  ;gc  iCv  nameis:  and  accress-.es,  irons;3. v:  vcni'Oring 

A  G  E i  C  :S~7R""  NL.MBER 


■;j.  distribution  availa3ii.7v  statement 


Approved  For  Public  Release;  Distribution  Unlimited 


Vx  r,,'r>4,lJ C  .vC'CS,' 

j  T*16  energy  spectrum  of  a  5  curie  Pu-Be  neutron  source  is  determined  bv  use  of  a 
i  neutron  time-of-f light  spectrometer.  This  spectrometer  measures  the  time  it  takes  for 
|  a  neutron  to  travel  between  two  plastic  scintillator  detectors.  Using  a  fixed  dis- 
j  tance  for  travel,  a  time  spectrum  is  generated  which  is  then  unfolded  into  the  energy 
i  spectrum  of  the  source.  The  different  aspects  of  equipment  selection,  connection 
!  snd  calibration  of  the  system  are  described.  With  the  resolved  spectrum  the  energv 
j  resolution  of  the  spectrometer  is  determined  by  using  a  Monte  Carlo  simulation  of  the 
1  experimental  set-up.  A  dual-parameter  experiment  is  conducted  in  which  the  proton 
i  recoil  spectrum  from  the  first  detector  is  placed  in  coincidence  with  the  time— of- 
j  flight  data.  Comparisons  with  the  published  spectrum  from  Arun  Kumar  and  P.  Nagar- 
!  ajan  are  made. 

rr 


ill  iUBDtCT  *2 3. MS 

4  Neutron  Spectrometers ,^Time^of“ Flight  Spectrometers,  Neutron 
energy  spectrum,  Pu-Be  Neutron  Source  Energy  spectrum,  <a _ 


• '  '...'.'.rT t  ::  -  c,\ <  • 

'  r  w  ? c< ~ 

_  lassified  I 

Aar 


-i  i-cjaiT  -  r<-.: 

■?  :  qs  :  •• 

Unclassified 


CN 


i£C-.rr '  i 

:f  as jt.-; Ar- 

Unclassified 


3  .UVBER  Cr 

]02 _ 

;  '.ON  or  A3,*  . 

SAR 


